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Future  high  speed  optical  fiber  communications  and  signal  processing  systems 
require  the  monolithic  integration  of  a laser  source  and  other  passive  as  well  as  active 
devices.  Key  to  integration  is  the  coupling  of  light  from  the  laser  source  to  a waveguide 
circuit  that  make  the  light  available  for  other  optical  processing  functions  such  as  modula- 
tion, switching  and  routing.  Existing  integration  schemes  either  require  complicated  and 
costly  material  regrowth  or  impose  severe  design  and  fabrication  restrictions.  This  disser- 
tation studies  a novel  monolithic  integration  of  a single  quantum  well  (SQW)  laser  and  a 
tapered  waveguide  coupler  in  the  GaAs/AlGaAs  material  system.  The  integrated  structure 
consists  of  two  sections:  the  laser  and  the  tapered  waveguide  coupler.  The  SQW 
waveguide  for  the  laser  and  the  multiple  quantum  well  (MQW)  waveguide  as  the  output 
waveguide  are  vertically  stacked.  In  the  laser  section,  the  MQWs  are  disordered  by  Zn 
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implantation  and  annealing  thus  becoming  the  cladding  layer  of  the  SQW  laser.  In  the 
tapered  waveguide  coupler  section,  the  MQWs  are  gradually  disordered  along  the 
waveguide  by  nitrogen  implantation  and  annealing.  This  permits  the  laser  output  couples 
to  the  MQW  waveguide  for  optical  processing  such  as  modulation.  A groove  etched 
between  the  laser  and  the  taper  sections  functions  as  a partially  transmissive  mirror  pro- 
viding the  feedback  for  the  laser.  The  tapered  waveguide  coupler  employs  an  adiabatic 
power  transfer  process  and  allows  for  independent  optimization  of  individual  components 
to  be  integrated.  The  vertical  configuration  of  the  two  waveguides  and  the  controlled 
impurity-induced  disordering  eliminate  the  regrowth  requirements  without  sacrificing  the 
device  design  and  fabrication  flexibility. 

In  this  dissertation,  the  key  technologies  for  the  integration  of  the  SQW  laser  and 
the  tapered  waveguide  coupler  have  been  investigated  theoretically  and  experimentally. 
The  MQWs  in  the  laser  section  are  disordered  by  Zn  implantation  and  subsequent  anneal- 
ing. A finite  difference  approach  has  been  developed  to  determine  the  concentration  (Zn) 
dependent  interdiffusion  and  the  model  leads  to  a prediction  of  a sharp  reduction  in  thresh- 
old current  density  due  to  the  disordering  of  the  MQWs  in  the  laser  section,  which  agreed 
well  with  our  experimental  results.  In  addition,  dependence  of  the  lasing  condition  due  to 
the  groove  are  investigated.  The  fabrication  techniques  including  the  MQWs  disordering 
by  implantation  and  annealing  and  the  reactive  ion  etching  (RIE)  for  the  groove  have  been 
developed.  Utilizing  these  techniques,  the  integration  of  an  SQW  laser  and  a tapered 
waveguide  coupler  in  the  GaAs/AlGaAs  material  system  has  been  accomplished. 
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CHAPTER  1 
INTRODUCTION 


The  breakthrough  in  semiconductor  diode  lasers  in  the  1960s  and  the  development 
of  low  loss  optical  fiber  in  the  1970s  were  responsible  for  the  explosion  in  the  fiber  optical 
communication  and  optical  signal  processing  systems.  The  progress  in  these  optical  sys- 
tems, in  turn,  has  stimulated  the  growth  of  a new  class  of  passive  and  active  optical  com- 
ponents in  thin  dielectric  layers  using  guided  wave  optics.  The  highly  promising  and 
sophisticated  technology  responsible  for  the  new  class  of  devices  is  known  as  integrated 
optics.  In  addition  to  the  advantages  of  greater  bandwidth  and  immunity  from  electromag- 
netic interference,  integrated  optic  devices  incorporate  increased  complex  circuitry  and 
functionality  along  with  high  reliability.  Successful  future  deployment  of  optical  fiber 
communications  and  optical  signal  processing  systems  relies  on  the  monolithic  integration 
of  a laser  source  with  various  passive  and  active  devices.  One  of  the  key  integrations  is  the 
coupling  of  light  from  the  laser  source  to  a waveguide  circuit  involving  other  optical  pro- 
cessing such  as  modulation,  switching  and  routing.  In  this  dissertation,  we  theoretically 
and  experimentally  investigate  a novel  laser-tapered  waveguide  coupler  integration;  in  par- 
ticular, we  develop  the  key  fabrication  techniques,  and  demonstrate  the  monolithic  integra- 
tion in  the  AlGaAs/GaAs  material  system.  The  first  section  of  this  introductory  chapter 
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describes  the  motivation  for  this  work.  In  section  1 .2,  we  compare  various  existing  laser- 
waveguide  integration  structures  and  introduce  the  novel  structure  we  proposed.  The  orga- 
nization of  the  remaining  chapters  of  this  dissertation  is  presented  in  section  1.3. 

1.1  Motivation 

Optical  fiber  offers  transmission  capacity  of  many  orders  of  magnitude  great  than 
that  of  conventional  transmission  lines.  The  overall  bandwidth  that  is  available  for  optical 
transmission  is  in  the  order  of  THz.  This  enormous  capacity  can  be  accessed  by  using 
techniques  such  as  wavelength  division  multiplex  (WDM);  it  allows  for  the  implementa- 
tion of  many  new  services  such  as  HDTV  distribution,  broadband  data  highways  and 
video  telephony.  In  a WDM  system,  the  overall  data  stream  is  divided  into  a set  of  sub- 
streams each  having  a different  wavelength.  The  optical  signals  are  multiplexed  by  an 
optical  power  combiner  and  then  demultiplexed  by  optical  filters  at  the  receiver  ends.  Fig- 
ure 1 shows  a multiwavelength  optical  broadcast  network.  The  network  based  on  the 
LAMBDANET  architecture  [1]  is  composed  of  a cluster  of  N communication  nodes.  Each 
node  transmits  its  information  at  a unique  wavelength  and  each  node  receives  all  the  infor- 
mation from  the  other  nodes.  The  nodes  are  optically  coupled  together  using  a transmis- 
sive optical  star  coupler  at  a hub  location  which  broadcasts  information  at  all  wavelengths 
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NODE  #2 


Figure  1.1  The  LAMBDANET  multiwavelength  architecture. 
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An  example  of  a local  video  distribution  network  using  the  LAMBDANET  archi- 
tecture is  shown  in  Fig.  2 [1].  All  the  transmitters  (lasers)  are  located  in  the  head  end  sta- 
tion and  all  other  end  stations  receive  the  video  broadcast.  At  least  256  video  channels, 
each  at  a data  rate  of  approximately  155  Mbit/sec  or  incorporating  time-multiplexed  chan- 
nels within  the  155  Mbit/sec  stream  for  even  more  video  channels,  are  optically  multi- 
plexed together  to  form  16  high  speed  data  channels  at  2.5  Gbit/sec,  which  are  transmitted 
to  and  demultiplexed  at  the  receiving  end  stations. 

Such  a multiwavelength  multiplex/demultiplex  network  demands  optoelectronic 
integration  for  superior  performance  and  high  reliability.  Integrated  transmitters,  receivers, 
and  optical  filters  are  continuously  being  developed.  Transmitters  generally  include  lasers, 
optical  amplifiers,  electronic  driving/modulation  circuits  and/or  integrated  external  optical 
modulators.  For  wavelength  multiplex  applications,  transmitter  arrays  are  often  required. 
Receiver  arrays  with  on-chip  wavelength  demultiplexing  are  desired  as  well.  Integrated 
receivers  employing  the  polarization  diversity  coherent  detection  scheme  incorporate  a 
variety  of  devices  including  a tunable  laser,  mode  converters,  mode  splitters,  3 dB  cou- 
plers, and  detectors.  Recently  a tunable  128  channel  optical  filter  has  been  demonstrated 
[2].  This  enables  transmission  of  thousands  of  video  channels.  Such  a highly  sophisticated 
video  distribution  network  can  be  realized  easily  with  optoelectronic  integration,  which 
improves  the  efficiency  and  simplifies  the  packaging  as  well  as  alignment  of  the  optical 
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One  of  the  key  problems  in  the  development  of  the  transmitter  is  the  modulation  of 
the  laser.  In  the  example  of  the  video  distribution  network  using  a LAMBDANET  archi- 
tecture (Fig.  1.2),  16  channels  of  155  Mbit/sec  video  signals  are  time-multiplexed  together 
to  form  a 2.5  Gbit/sec  high  speed  data  stream  that  modulates  the  laser  to  generate  the  opti- 
cal signals  for  transmission.  High  degree  of  time-multiplexing  or  other  techniques  for 
increasing  the  number  of  channels  would  result  in  an  even  higher  speed  modulation  of  the 
laser.  However,  high  speed  direct  modulation  of  a semiconductor  diode  laser  causes  two 
problems:  dynamic  wavelength  shift  (also  known  as  frequency  chirping)  and  dynamic 
spectrum  broadening. 

Direct  modulation  of  a diode  laser  is  performed  by  modulating  the  injection  cur- 
rent. The  modulation  of  the  injection  current  causes  the  modulation  of  the  carrier  density 
in  the  active  layer  of  the  diode  laser.  It  is  well  known  that  the  change  in  the  carrier  density 
results  in  the  change  in  the  refractive  index  of  the  active  layer  due  to  the  change  of  the 
band-to-band  absorption.  A refractive  index  change  in  the  active  layer  means  a change  of 
the  effective  optical  length  in  the  laser  cavity,  and  therefore  results  in  a shift  of  the  lasing 
wavelength.  This  phenomenon  is  called  dynamic  wavelength  shift,  also  known  as  fre- 
quency chirping  [3-8].  On  the  other  hand,  as  the  carrier  density  in  the  active  layer  oscil- 
lates under  the  modulation,  accordingly  the  medium  gain  also  oscillates  in  concert.  During 
this  process,  the  mode  gains  of  various  longitudinal  modes  exceed  the  threshold  tran- 
siently so  that  multimode  lasing  occurs  causing  dynamic  spectral  broadening  [9-13]. 
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Dynamic  wavelength  shift  and  spectral  broadening  are  significant  under  high 
speed  direct  modulation.  The  wavelength  division  multiplexing  and  long  distance  optical 
fiber  transmission  require  single  wavelength  or  narrow  line  width  lasers. Dynamic  wave- 
length shift  and  spectral  broadening  will  reduce  the  bandwidth  of  the  WDM  and  optical 
transmission. 

To  avoid  these  shortcomings,  an  alternative  to  the  direct  modulation  of  the  diode 
laser  is  to  use  an  external  optical  modulator.  An  additional  advantage  of  using  an  external 
modulator  is  that  modulators  using  quantum  well  effects  may  have  higher  modulation 
speed  than  direct  modulation  of  the  laser  diode  whose  speed  is  limited  by  the  carrier  life 
time  during  the  recombination. 

Monolithic  integration  of  a laser  and  a modulator  has  been  an  active  area  of 
research  in  recent  years.  Double  heterostructure  or  quantum  well  lasers  and  multiple  quan- 
tum well  waveguide  modulators  have  been  used  in  various  integration  structures  [14].  A 
key  problem  in  the  integration  is  the  coupling  of  the  light  from  the  laser  to  a waveguide  to 
make  the  light  available  for  modulation.  Two  factors  should  be  considered  in  the  design  of 
the  integrated  structure.  To  make  the  monolithic  integration  practical  and  cost-effective, 
the  fabrication  process  must  be  as  simple  as  possible.  Among  semiconductor  processing 
technologies  available  today,  MOCVD/MBE  growth  is  among  the  most  sophisticated 
techniques  and  is  quite  expensive.  Therefore,  if  possible,  regrowth  should  be  avoided  in 
the  integrated  structure  to  minimize  complexity  in  fabrication.  The  optimization  of  the 
performance  of  the  integrated  structure  is  another  important  issue  in  the  design  of  the 
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Structure.  It  should  be  noted  in  this  connection  that  devices  to  be  integrated  on  the  same 
chip  may  have  conflicting  requirements  for  layer  compositions  and  dimensions.  In  view  of 
this,  it  would  be  desirable  to  have  an  integrated  structure  that  allows  for  independent  opti- 
mization of  individual  components. 

This  dissertation  concentrates  on  the  study  of  the  first  part  of  a laser-modulator 
transmitter  chip,  namely,  a novel  monolithic  integration  of  a single  quantum  well  (SQW) 
laser  and  a tapered  waveguide  coupler.  In  addition  to  avoiding  the  complicated  and  costly 
material  regrowth,  the  integration  scheme  of  a SQW  laser  and  a tapered  coupler  offers 
flexibility  for  design  and  optimization  of  individual  components  to  be  integrated.  Its  appli- 
cation is  not  limited  to  the  laser-modulator  integration  as  it  also  makes  the  light  output 
from  the  laser  available  through  a waveguide  coupler  for  other  optical  signal  processing 
applications  such  as  switching  and  routing.  In  the  following  section,  we  compare  various 
existing  integration  structures  and  introduce  the  novel  structure. 

1.2  The  Laser- Waveguide  Integrated  Structures 

For  the  laser-waveguide  integration,  feedback  for  the  lasing  action  and  coupling  of 
light  from  the  laser  to  the  waveguide  must  be  established  in  the  integrated  structure.  Con- 
ventional semiconductor  lasers  use  cleaved  crystal  facets  as  reflectors.  Cleaving  is  a sim- 
ple procedure  and  the  cleaved  facet  is  flat  to  an  atomic  scale.  However,  it  is  difficult  to 
apply  this  technique  in  a monolithic  integrated  laser  structure.  Etched  facets  and  distrib- 
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uted  Bragg  reflectors  (DBR)  are  usually  used  for  integrated  lasers.  Reflectors  using  etched 
facets  can  be  fabricated  by  wet  chemical  etching  or  dry  etching  such  as  reactive  ion  etch- 
ing (RIE).  Since  the  etched  facets  are  defined  by  the  photolithographic  processes,  they  can 
be  made  virtually  in  any  shape,  size,  and  crystal  orientation  and  thus  a variety  of  cavities 
can  be  fabricated.  A distributed  Bragg  reflector  in  a semiconductor  laser  usually  has  peri- 
odic corrugation  on  the  active  or  passive  optical  waveguide  layer.The  equivalent  refractive 
index  of  the  optical  waveguide  is  hence  modulated  periodically  along  the  waveguide  prop- 
agation direction  so  that  forward  propagating  light  is  scattered  backward  by  Bragg  reflec- 
tion if  the  wavelength  of  the  light  is  close  to  the  Bragg  wavelength.  The  periodic 
corrugations  for  the  distributed  Bragg  reflector  in  the  semiconductor  are  generally  first 
patterned  by  holography  or  e-beam  lithography  and  then  etched  by  either  wet  chemical  or 
dry  etching.  A laser  with  periodic  corrugations  on  the  passive  optical  waveguides  at  two 
ends  of  the  active  waveguide  is  referred  to  as  a distributed  Bragg  reflector  (DBR)  laser.  A 
DBR  laser  is  wavelength  tunable  by  the  modulation  of  the  refractive  index  of  the  distrib- 
uted Bragg  reflector.  On  the  other  hand,  the  periodic  corrugation  can  be  fabricated  directly 
on  the  active  optical  waveguide  to  form  a distributed  feedback  laser  (DFB).  With  the 
advantage  of  simplicity  over  DBR  lasers,  DFB  lasers  are  commonly  used  for  integration. 

To  date,  two  types  of  structures  for  coupling  light  from  the  laser  to  the  waveguide 
have  received  widespread  use.  One  is  butt-coupling,  and  the  other,  directional/evanescent 
coupling.  In  a butt-coupling  structure,  the  active  waveguide  is  directly  coupled  to  the  out- 
put passive  waveguide.  Usually  a DFB  laser  is  used,  in  which  an  etched  grating  on  the 
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Figure  1.3  DFBlaser  -waveguide  integration  using  the  butt-coupling 
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active  waveguide  functioning  as  a distributed  Bragg  reflector  provides  reflection  for  the 
lasing  action.  Figure  1.3  shows  a typical  butt-coupling  scheme  [15,  16].  Three  growths  are 
required  to  fabricate  such  a structure.  The  lower  cladding  layer  and  the  active  waveguide 
layers  are  first  grown  on  the  substrate.  After  the  grating  is  etched  over  the  active 
waveguide  using  holography/e-beam  lithography  and  wet/dry  etching,  a regrowth  is  per- 
formed to  form  the  upper  cladding  and  ohmic  contact  layers.  Generally,  the  output  passive 
waveguide  has  material  compositions  and  dimensions  different  from  those  of  the  laser 
waveguide.  To  accommodate  this  requirement,  the  active  waveguide  outside  the  laser  por- 
tion is  etched  away  and  another  regrowth  needs  to  be  performed  to  form  the  passive  output 
waveguide.  The  main  problem  with  this  approach  is  the  cost-effectiveness  and  the  com- 
plexity of  fabrication.  Additionally,  the  axial  offset  between  the  active  and  passive 
waveguides  must  be  minimized  by  accurate  control  of  the  etching  and  regrowth  process. 
When  the  optimizations  of  the  laser  waveguide  and  the  component  to  be  integrated  such  as 
a modulator  lead  to  different  dimensions  in  the  active  and  passive  waveguides,  mode 
matching  and  hence  the  coupling  efficiency  between  the  two  waveguides  may  become  a 
significant  problem. 

In  the  directional  coupling  scheme  for  the  laser-waveguide  integration,  the  active 
and  passive  waveguides  are  placed  in  close  proximity  so  that  the  mode  fields  of  the  indi- 
vidual waveguides  overlap.  The  working  principle  is  the  same  as  that  of  a directional  cou- 
pler. When  light  is  launched  into  one  of  the  two  identical  waveguides  of  a directional 
coupler,  generally  a symmetric  (fundamental)  normal  mode  and  an  antisymmetric  (first 
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Active  Waveguide 


Output  Waveguide 


Substrate 


(b) 


Figure  1.4  (a)  Directional  coupler  and  (b)  Integrated  twin  guide  (ITG) 
laser  using  directional  coupling. 
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order)  normal  mode  are  excited.  The  two  normal  modes  are  in  phase  at  the  input  end  and 
the  fields  of  the  two  normal  modes  add  constructively  in  the  guide  A and  destructively  in 
the  guide  B so  the  light  is  confined  in  the  input  waveguide  A.  The  two  normal  modes  have 
different  propagation  constants  Pq  and  Pj,  respectively.  After  propagating  a distance  of  the 
coupling  length,  the  two  modes  are  1 80°  out  of  phase  and  add  destructively  in  the  guide  A 
and  constructively  in  the  guide  B so  the  power  is  confined  in  the  guide  B.  Essentially,  the 
power  is  transferred  from  the  guide  A to  the  guide  B over  a distance  of  the  coupling  length 
/=7c/(Po-Pi).  The  two  individual  waveguides  may  have  different  refractive  indices  and 
dimensions  but  they  must  be  designed  so  as  to  have  identical  propagation  constants  at  the 
operating  wavelength  to  achieve  complete  power  transfer  between  the  two  guides.  Also, 
the  length  of  the  coupling  region  should  be  equal  to  the  coupling  length  which  is  deter- 
mined by  material  compositions  and  waveguide  dimensions. 

The  integrated  twin  guide  (ITG)  laser  structure  shown  in  Figure  1.4  employs  the 
directional  coupling  scheme  [17-20].  The  laser  active  waveguide  and  the  passive  output 
waveguide  are  configured  as  a vertical  directional  coupler  with  light  being  coupled  from 
the  laser  waveguide  to  the  output  waveguide  below.  Etched  facets  in  the  active  waveguide 
are  used  as  reflectors  for  the  laser.  The  structure  can  be  fabricated  in  a single  growth.  How- 
ever, due  to  the  interferometric  nature  of  the  directional  coupler  as  we  discussed  above,  the 
two  waveguides  must  have  identical  propagation  constants  to  achieve  efficient  coupling. 
This  imposes  severe  restrictions  on  the  device  design  and  strict  tolerance  on  the  fabrica- 
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tion.  Integrated  laser  structures  using  directional/evanescent  coupling  tend  to  have  higher 
threshold  current  since  the  optical  field  is  distributed  between  the  two  waveguides  result- 
ing in  a lower  optical  confinement  factor  in  the  active  layer. 

To  overcome  the  shortcomings  of  these  structures,  recently,  a novel  monolithic 
integration  of  an  SQW  laser  and  an  MQW  modulator  via  a tapered  waveguide  interconnect 
in  the  AlGaAs/GaAs  material  system  (Figure  1.5)  has  been  proposed  and  studied  in  our 
group  [21-23].  The  integrated  structure  consists  of  three  sections:  laser,  tapered  coupler, 
and  modulator.  The  key  technology  is  the  integration  of  the  SQW  laser  and  the  tapered 
waveguide  coupler  which  couples  the  light  from  the  SQW  laser  to  a MQW  waveguide  for 
modulation.  In  this  integrated  structure,  all  layers  can  be  grown  in  a single  growth.  Unlike 
the  uniform  directional  coupler,  the  tapered  coupler  employs  an  adiabatic  power  transfer 
process  and  does  not  require  the  two  waveguides  to  have  identical  propagation  constants 
over  the  entire  length.  Somewhat  relaxed  design  constraints  and  fabrication  tolerance  as 
well  as  the  ability  to  separately  optimize  individual  components  provide  definite  advan- 
tages to  this  approach.  In  this  dissertation,  we  focus  on  the  study  of  the  monolithic  integra- 
tion of  the  SQW  laser  and  the  tapered  waveguide  coupler.  The  integrated  structure  and  its 
working  principles  will  be  discussed  in  detail  in  Chapter  2 and  key  technologies  will  be 
further  studied  in  the  remaining  chapters.  The  organization  of  this  dissertation  is  presented 


in  the  next  section. 
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Figure  1 .5  A novel  monolithic  integration  of  an  SQW  laser  and  an  MQW 
modulator  via  a tapered  waveguide  coupler. 
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1.3  Organization  of  the  Dissertation 

Following  the  introduction  presented  in  this  chapter,  Chapter  2 describes  the  struc- 
ture and  working  principle  of  the  novel  monolithic  integration  of  a SQW  laser  and  a taper 
waveguide  coupler.  We  start  with  the  discussion  of  multi-layer  optical  waveguides  and 
quantum  well  lasers  which  form  the  basic  components  of  the  integrated  structure.  Then  the 
principle  of  the  tapered  waveguide  coupler  comprising  of  a single  quantum  well  (SQW) 
waveguide  and  a multiple  quantum  well  (MQW)  waveguide  will  be  illustrated.  The  last 
section  of  Chapter  2 discusses  the  principle  of  operation,  fabrication  process,  and  the 
advantages  of  the  entire  integrated  structure. 

The  key  technologies  for  the  fabrication  of  the  novel  integrated  structure  are  fur- 
ther studied  in  the  remaining  chapters.  Three  areas  of  research  including  the  disordering  of 
the  MQWs,  etched  groove,  and  the  fabrication  technologies  are  covered.  To  achieve  low 
threshold  current  for  the  SQW  laser,  it  is  necessary  to  disorder  the  MQWs  in  the  laser  sec- 
tion. In  Chapter  3,  we  numerically  model  the  MQW  disordering  process  by  Zn  and  its 
influences  on  the  threshold  current  of  the  laser.  The  dependences  of  the  Zn  profile,  the 
compositional  and  refractive  index  profiles  of  the  disordered  MQWs,  the  optical  confine- 
ment factor  and  the  laser  threshold  current  on  the  disordering  parameters  are  calculated. 
Chapter  4 presents  the  relations  between  the  reflectance  of  the  groove,  the  laser  threshold 
current,  the  output  power,  and  the  groove  parameters.  In  Chapter  5,  we  study  the  key  fabri- 
cation technologies  including  the  MQWs  disordering  by  Zn  implantation  and  annealing 
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and  the  groove  etching  by  reactive  ion  etching  (RIE).  Experimental  results  are  compared 
with  those  of  the  numerical  modeling  and  the  integration  is  demonstrated.  In  the  last  chap- 
ter, summary  and  conclusions  of  this  dissertation  are  presented. 


CHAPTER  2 

NOVEL  INTEGRATION  OF  A SQW  LASER  AND  A TAPERED  COUPLER 


This  chapter  describes  and  analyzes  the  novel  monolithic  integration  of  a single 
quantum  well  (SQW)  laser  and  a tapered  waveguide  coupler.  We  first  present  the  basic 
optical  components  used  in  this  integration,  namely,  multiple  layer  optical  waveguide, 
quantum  well  laser,  and  tapered  waveguide  coupler.  Then  we  describe  the  integrated  struc- 
ture, its  principle  of  operation,  and  the  techniques  employed  in  fabrication  of  the  structure 
as  well  as  its  advantages  over  other  integrated  structures. 

2.1  Multilayer  Optical  Waveguide 

The  optical  waveguides  are  the  building  blocks  commonly  used  in  both  active  and 
passive  integrated  optic  devices  where  the  light  is  confined  and  guided  in  the  waveguide.  A 
well  known  example  of  the  optical  waveguides  is  the  optical  fiber  which  has  a higher 
refractive  index  core  surrounded  by  a low  refractive  index  cladding  with  usually  a circular 
cross  section.  Waveguides  in  integrated  optic  devices  can  have  planar,  channel,  ridge  and 
rib  structures.  Figure  2. 1 shows  a simple  slab  planar  waveguide  which  consists  of  three 
layers,  namely,  a substrate,  a film,  and  a cover  with  refractive  indices  of  n^,  nf,  and  n^,. 
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Figure  2. 1 A three  layer  slab  waveguide. 
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respectively.  To  guide  light  in  the  film,  nf  must  be  greater  than  n^  and  n^..  The  cover  layer 
may  be  air  in  some  devices.  The  layered  structure  can  be  formed  by  material  deposition, 
diffusion,  or  epitaxial  growth. 

A ray  optical  model  of  light  propagation  in  this  simple  waveguide  can  be  used  to 
illustrate  the  basic  concepts  of  optical  waveguides.  The  underlying  principle  of  guidance 
in  the  optical  waveguide  is  the  total  internal  reflection  at  the  interfaces.  When  light  is  inci- 
dent on  an  interface  between  two  media,  reflection  and  refraction  occur.  The  reflectance 
and  the  refractive  angles  obey  Snell’s  and  Fresnel’s  laws,  respectively.  If  the  light  is  inci- 
dent on  the  interface  from  the  higher  refractive  index  (ny)  medium  to  the  lower  index  (n^  or 

n^)  medium  and  the  incident  angle  is  greater  than  the  so-called  critical  angle  sin^(nj/ng  ^), 
total  internal  reflection  occurs.  The  amplitude  of  the  reflectance  equals  one  and  a phase 
shift  2<^  is  induced  at  the  interface.  The  reflectance  can  be  expressed  as 

r = expO'2(j)) 

The  phase  shift  2(])  due  to  total  internal  reflection  is  determined  by  the  incident  angle  0 and 
the  refractive  indices. 

The  three  layer  waveguide  shown  in  Figure  2.1  has  two  critical  angles,  for  total 
reflection  at  the  film-substrate  interface  and  6^  (<65)  at  the  film-cover  interface.  We 
assume  nf>ns>nj..  When  the  incident  angle  0 is  smaller  than  0^  and  0^.,  total  reflection  does 
not  occur  at  the  two  interfaces  and  the  light  incident  from  the  substrate  is  refracted  through 
the  film  and  cover  layers  and  escapes  from  the  guide  (Figure  2.2a).  There  is  no  confine- 
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c.  Guided  mode  at  discreet  angles 


Figure  2.2  Modes  in  the  slab  waveguide 
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ment  and  guiding  of  the  light.  This  mode  of  light  propagation  is  called  a true  radiation 
mode.  When  the  incident  angle  is  increased  such  that  0£.<0<0s,  the  total  reflection  occurs 
at  the  film-cover  interface  but  refracted  at  the  film-substrate,  as  shown  in  Figure  2.2b.  The 
light  again  escapes  from  the  guide  on  the  substrate  side.  This  is  called  a substrate  radiation 
mode.  When  the  incident  angle  is  further  increased  to  0>0j>0j.,  total  reflection  occurs  at 
both  interfaces.  Under  this  condition,  light  entering  the  film  layer  may  be  confined  and 
guided  in  a zig-zag  path  in  the  film  (Figure  2.2c)  and  we  have  a guided  mode. 

However,  for  the  guided  mode,  not  all  angles  of  incidence  0 are  allowed  even 
though  they  are  greater  than  both  critical  angles  0^  and  0(..  When  the  light  propagates  in  a 
zig-zag  path  along  the  waveguide,  the  total  transverse  phase  shift  (1)^  must  satisfy  the  self- 
consistency  condition  (also  called  the  transverse  resonance  condition),  which  is  (1)t=2v7U 
(v=0,l,2...).  The  phase  shifts  in  this  three  layer  waveguide  include  the  shift  due  to  a round- 
trip  transverse  passage  through  the  film,  2knfhcos0,  and  shifts  due  to  the  total  internal 
reflections  at  the  two  interfaces,  -2(t)s  and  -2^c-  The  self-consistency  condition  becomes 

2^«y/jcos0-2(t)^-2(l)^  = 2v7C 

where  v=0,  1,  2,...  Note  that  (t)^  and  (])(.  are  also  functions  of  0.  This  condition  shows  that 
only  a discrete  set  of  angles  0 lead  to  guided  modes.  The  angles  with  v=0,  1,  2,...  corre- 
spond to  the  so-called  fundamental  mode,  first-order  mode,  second-order  mode,...,  respec- 
tively. From  the  above  equation,  it  is  possible  to  make  the  waveguide  single-moded  (the 
equation  holds  only  with  v=0)  by  properly  choosing  the  values  of  h,  nj,  Uj,  and  n^,.  Once  0^ 
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is  determined  from  the  self-consistency  condition,  the  propagation  constants  of  the  guided 
modes  can  be  calculated  from  ^^=knjsin%.  Thus,  the  self-consistency  condition  also  gives 
the  dispersion  relation  of  the  waveguide.  This  dispersion  relation  as  well  as  the  mode  field 
profiles  can  also  be  obtained  by  using  the  electromagnetic  theory  of  waveguides. 

In  practical  applications,  optical  waveguides  often  have  multilayered  structures. 
Semiconductor  separate  confinement  heterostructure  (SCH)  lasers  use  a five  layer 
waveguide  to  achieve  separate  confinement  for  the  light  and  the  injected  carriers[24].  The 
graded  index  (GRIN)  SCH  structure  used  for  single  quantum  well  (SQW)  lasers  to  maxi- 
mize optical  confinement  can  be  treated  numerically  as  a multilayer  waveguide[25,  26]. 
Many  guided  wave  devices  that  consist  of  two  or  more  single  mode  waveguides  such  as 
vertical  directional  couplers  can  also  be  considered  as  a multilayer  waveguide.  In  fact,  the 
two  components  in  our  integrated  structure,  the  SQW  laser  and  the  tapered  coupler,  are 
both  constructed  from  the  multilayer  waveguide  structure.  The  basic  principles  for  the 
confinement  and  guiding  of  light  in  a multilayer  waveguide  are  the  same  as  those  in  a sim- 
ple three  layer  waveguide.  However,  more  complicated  techniques  such  as  the  multilayer 
stack  theory[27]  are  used  to  analyze  the  dispersion  relation  and  field  profiles  of  the  guided 
modes  in  a multilayer  waveguide.  An  analysis  of  the  tapered  coupler  as  a multilayer 
waveguide  is  presented  in  section  2.3  and  the  multilayer  stack  theory  is  employed  in  sec- 
tion 3.4  to  model  the  optical  confinement  dependence  of  the  SQW  laser  on  the  MQW  dis- 
ordering. 
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2.2  Quantum  Well  Lasers 

In  this  section,  we  review  the  basics  of  the  quantum  well  laser  which  is  one  of  the 
key  components  in  our  integration.  Early  semiconductor  lasers  were  just  a simple  GaAs  p- 
n homojunction  in  a Fabry-Perot  cavity  formed  by  two  cleaved  facets  where  injected  hole 
and  electron  carriers  recombine  and  emit  photons[28,  29].  The  lasing  thresholds  of  these 
lasers  were  high  and  they  could  operate  only  in  pulsed  mode.  The  development  of 
AlGaAs/GaAs  double  heterostructure  (DH)  lasers  made  the  room  temperature  CW  opera- 
tion possible[30-32].  The  bandgap  difference  between  AlGaAs  and  GaAs  provides  barri- 
ers in  the  DH  structure  that  confine  injected  carriers,  thus  increasing  the  carrier  density. 
The  refractive  index  difference  between  the  two  materials  makes  the  DH  structure  an  opti- 
cal waveguide  that  confines  the  photons  increasing  the  overlap  of  the  optical  field  and  the 
gain  medium.  The  confinements  of  the  injected  carriers  and  the  photons  greatly  reduce  the 
lasing  threshold.  The  separate  confinement  heterostructure  (SCH)  optimizes  the  confine- 
ments of  the  carriers  and  photons[33],  helps  in  placing  the  corrugated  grating  away  from 
the  active  layer  in  the  DFB  lasers[34]  and  is  a necessity  for  quantum  well  lasers  in  order  to 
maximize  the  optical  confinement. 

With  their  superior  in  performance,  quantum  well  lasers  are  now  replacing  the  con- 
ventional DH  lasers.  We  begin  the  comparison  of  the  quantum  well  and  the  conventional 
bulk  lasers  from  the  density  of  states.  In  a semiconductor  crystal,  the  electronic  states  in 
the  conduction  and  valence  bands  are  described  by  the  electron  wavefunctions  which  are 
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governed  by  the  Schroedinger  equation 


+V(r) 

2m 


V|/  = E\\f 


where  V(r)  is  the  potential  created  by  the  crystal  lattice,  \ji  is  the  wavefunction  and  E the 
electronic  energy  associated  with  the  wavefunction.  The  wavefunctions  are  uniquely  rep- 


resented by  Bloch  wave  \|/  = e^’^u{k,r)  and  characterized  by  the  electron  wavevector  k. 
The  fact  that  the  electron  wavefunction  must  go  to  zero  at  the  edges  of  the  bulk  crystal 
gives  the  boundary  conditions  for  the  wavefunctions  resulting  in  a discrete  set  of  k vectors 
k = n^kj^  + Hyky  + n^k^  where  rii  = 0,  ±1,  ±2,...,  |/:,|  = n/L^  and  L,-  are  the  crystal 

dimensions  with  i=x,  y,  and  z-  The  wavefunctions  with  these  k vectors  represent  the  possi- 
ble states  an  electron  can  have.  In  III-V  semiconductors,  the  energy  E and  the  k vector  of 

h^k^ 

an  electron  are  generally  related  by  £ = E^  + — — , the  so-called  parabolic  band  structure 

approximation.  Figure  2.3.  la-c  show  a cubic  bulk  crystal,  the  k vector  states  in  the  k space 
and  the  energy-^  vector  relation,  respectively.  Note  that  each  cross  point  in  the  k vector 
space  represents  a possible  electronic  state  and  the  distance  between  two  adjacent  points  is 
|kj|  = n/L-,  which  is  dependent  on  the  crystal  dimensions.  If  we  reduce  one  of  the 


dimensions  («L^  and  Ly),  the  distance  between  the  two  adjacent  points  in  the  z direction 
is  enlarged  and  the  uniformly  distributed  k states  become  far-separated  planes  of  states 
corresponding  to  different  quantum  numbers  n^.  Accordingly,  the  energy  is  also  divided 
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Figure  2.3  The  k state  space,  E-k,  and  p(E)  for  bulk  crystal  and  quantum  well. 
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into  subbands.  These  are  illustrated  in  Figure  2.3.2a-c.  When  is  less  than  1000  A for  III- 
V semiconductors,  the  quantum  effect  becomes  prominent  and  the  structure  is  called  a 
quantum  well.  For  a quantum  well,  is  many  orders  of  magnitude  larger  than  and  ky 
The  problem  of  the  density  of  states  becomes  two  dimensional  in  each  plane  and  the 
planes  of  states  in  the  k space  are  far  apart. 

With  the  help  of  Figures  2.3.1b  and  2.3.2b,  the  density  of  states  in  the  k space  can 
be  calculated  for  the  bulk  and  quantum  well  structures.  From  the  relation  between  the 
energy  E and  the  k vector,  the  density  of  states  in  term  of  energy  can  be  obtained.  After 
extensive  and  tedious  derivations,  it  can  be  shown[35]  that  for  the  bulk  crystal,  the  density 
of  states  p(£)  is  parabolic-like  while  for  the  quantum  well,  p(£)  is  staircase-like  discontin- 
uous function  with  constant  density  of  states  at  each  subband,  as  shown  in  Figures  2.3.  Id 
and  2.3.2d,  respectively.  It  should  be  noted  that  the  staircase-like  density  of  states  of  a 
quantum  well  will  reduce  to  that  of  the  bulk  crystal  with  increasing  well  thickness. 

It  is  the  staircase-like  density  of  states  of  the  quantum  well  that  makes  the  perfor- 
mances of  the  quantum  well  lasers  superior  to  those  of  the  conventional  bulk  active  layer 
DH  laser.  Photon  emission  in  a semiconductor  laser  is  the  result  of  the  recombination  of 
the  injected  electrons  and  hole  carriers.  The  carrier  density  n(E)  is  the  product  of  the  den- 
sity of  states  p(£)  and  the  Fermi-Dirac  distribution  function  f{E)  which  determines  the 
probability  of  an  carrier  occupying  an  energy  level  E and  is  expressed  as 

1 

l + exp[(£-Ep/(kgT)] 


fm  = 
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where  Ef  is  the  quasi-Fermi  level  determined  by  the  total  number  of  injected  carriers.  The 
p(£),  j{E)  and  n{E)  for  the  bulk  crystal  and  the  quantum  well  are  plotted  in  Figures  2.4a 
and  2.4b,  respectively.  The  differences  in  the  carrier  density  profiles  as  a result  of  the  dif- 
ferent densities  of  states  lead  to  different  responses  to  the  increase  in  the  total  number  of 
injected  carriers.  In  the  bulk  crystal,  the  carrier  density  as  a function  of  energy  increases 
gradually  since  the  density  of  states  increases  with  increasing  energy.  It  reaches  a maxi- 
mum and  then  decreases  because  the  electrons  have  lower  probability  of  occupying  higher 
energy  levels  as  the  Fermi-Dirac  distribution  decreases  with  increasing  energy.  For  the 
quantum  well,  because  of  the  constant  value  of  the  density  of  states,  electrons  occupy  the 
bandedge  immediately  with  the  peak  of  the  carrier  density  located  at  the  bandedge  with  a 
narrower  width.  This  leads  to  a much  larger  increase  in  the  bandedge  carrier  density  in 
response  to  the  increase  of  the  total  number  of  injected  carriers  which  corresponds  to  a 
shift  of  the  Fermi  level  to  a higher  energy  level.  The  enhancement  of  the  bandedge  carrier 
density  greatly  increases  the  peak  gain  making  a low  lasing  threshold  possible.  It  should 
be  noted  that  with  a quantum  wire  structure,  the  density  of  states  decreases  with  increasing 
energy  level  and  the  bandedge  carrier  density  will  be  further  enhanced. 

With  the  advances  in  the  epitaxial  growth  technology,  the  quantum  well  lasers 
were  demonstrated  experimentally  and  their  advantages  exemplified  [36-38].  For  the 
GaAs/Al^Caj.^As  material  system,  AlGaAs  that  has  a larger  bandgap  is  used  as  the  barrier 
while  a lower  bandgap  GaAs  or  low  A1  concentration  (small  x)  AlGaAs  layer,  sandwiched 
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a.  bulk 


b.  quantum  well 


Figure  2.4  schemetic  plots  of  n(E),  r(E)  and  f(E). 
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between  two  barrier  layers,  acts  as  the  quantum  well.  Typically,  the  well  thickness  is 
below  100  A.  Such  a thin  layer  of  high  quality  can  be  grown  either  by  molecular  beam  epi- 
taxy (MBE)  or  by  metal  organic  chemical  vapor  deposition  (MOCVD).  For  quantum  well 
lasers,  it  is  necessary  to  use  a separate  confinement  heterostructure  (SCH)  to  improve  the 
optical  confinement  in  the  active  well  layer.  Step  index  or  graded  index  (GRIN)  SCH  is 
employed  to  form  an  optical  waveguide  in  order  to  confine  the  lasing  light  [39].  In  Section 
3.5,  we  will  investigate  the  dependence  of  the  SQW  laser  threshold  current  on  the  MQW 
disordering  parameters. 


2.3  Tapered  Waveguide  Coupler 

Tapered  directional  coupler  is  another  key  component  in  the  laser-waveguide  inte- 
gration. It  couples  the  SQW  laser  waveguide  to  the  MQW  output  waveguide.  Although  the 
structure  of  a tapered  coupler  is  similar  to  that  of  a conventional  directional  coupler,  its 
principle  of  operation  is  completely  different.  Figure  2.5  shows  a tapered  directional  cou- 
pler that  consists  of  two  planar  waveguides  with  the  top  guide  thickness  tapered.  Curves 
labeled  a and  b in  Figure  2.5  represent  the  propagation  constants  of  the  two  individual 
waveguides.  The  two  modes  are  synchronous  at  the  center  of  the  coupler  where  the  two 
waveguides  have  identical  thicknesses.  The  power  transfer  between  the  two  waveguides 
can  be  analyzed  using  the  concept  of  local  normal  modes.  We  mentioned  in  Section  2. 1 
that  the  structure  of  a two-guide  coupler  can  be  viewed  as  a multilayer  waveguide.  In  this 
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Figure  2.5  (a)  A tapered  directional  coupler  with  top  guide 
thickness  tapered,  (b)  Propagation  constants  of  the  modes  of 
individual  guides  and  (3^  and  local  normal  modes  Pq  and  Pj 
of  the  composite  structure. 
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case,  the  composite  structure  is  a five  layer  waveguide.  Since  the  waveguide  is  tapered,  the 
mode  characteristics  will  vary  along  the  five-layer  waveguide  and  therefore,  we  must  deal 
with  the  local  normal  modes.  The  local  normal  modes  are  the  normal  modes  of  the  five 
layer  waveguide  at  a given  position  along  its  length.  Using  the  multilayer  stack  theory 
[27],  the  propagation  constants  Pq  ^^d  Pi  and  corresponding  field  profiles  of  the  local  nor- 
mal modes  as  a function  of  the  position  can  be  calculated.  The  Po  and  Pi  are  illustrated  in 
Figure  2.5  and  the  field  profiles  as  a function  of  the  positions  along  the  length  of  the 
slowly  varying  taper  are  also  shown  in  the  figure.  Note  that  at  the  input  end  of  the  coupler, 
the  fundamental  mode  (Pq)  has  its  power  concentrated  in  the  upper  guide  and  at  the  output 
end,  the  power  is  concentrated  in  the  lower  guide.  For  the  first  order  mode  (Pi),  the  power 
shifts  from  the  lower  guide  to  the  upper  guide.  Essentially,  due  to  the  gradual  tapering,  the 
normal  modes  evolve  slowly  along  the  coupler  altering  the  field  profiles  and  resulting  in 
power  transfer  between  the  two  waveguides.  If  we  excite  only  one  normal  mode  in  the 
coupler,  we  would  expect  to  achieve  a complete  power  transfer  from  one  guide  to  the 
other.  However,  an  important  condition  must  be  met.  The  mode  propagation  along  the 
tapered  coupler  can  be  analyzed  using  the  step  transition  model  [40, 41]  in  which  the  con- 
tinuously tapered  structure  is  approximated  by  a series  of  discontinuous,  abrupt  steps. 
Since  the  local  normal  modes  in  two  adjacent  steps  are  not  orthogonal,  coupling  between 
them  can  occur  at  each  step  meaning  the  power  in  the  fundamental  mode  can  be  coupled 
into  the  first  order  mode  in  the  next  step.  This  would  cause  large  crosstalk  in  the  coupler. 
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However,  if  the  tapering  of  the  waveguide  is  sufficiently  gradual,  the  coupling  between  the 
two  local  normal  modes  is  negligible  and  the  power  will  remain  in  one  mode  and  trans- 
form itself  from  being  concentrated  in  one  guide  to  being  concentrated  in  the  other.  This  is 
the  adiabatic  operation.  An  analysis  [41]  using  the  step  transition  model  yields  three  con- 
ditions for  the  tapered  coupler: 

1)  < 1.5  for  adiabatic  operation,  (2.1) 

KK* 

2)  AZj-  > 213/(Apj-)  for  complete  power  transfer  (>95%),  (2.2) 

and 

3)  large  (AP)/|k|  ensuring  only  one  local  normal  mode  is  excited  by  power 

launched  into  one  guide,  (2.3) 

where  AP^-  is  the  total  variation  of  the  difference  in  the  propagation  constants  of  the 

tapered  guide  over  the  taper  length  AZj- , and  K is  the  coupling  coefficient  between  the  two 
individual  guides. 

These  conditions  show  the  advantages  of  the  tapered  couplers  over  the  conven- 
tional directional  couplers.  Recall  our  brief  discussion  of  the  latter  in  Section  1.2.  The  nor- 
mal mode  field  profiles  do  not  change  when  the  modes  propagate  along  the  conventional 
coupler  because  it  uses  two  uniform  waveguides.  The  power  transfer  is  due  to  the  con- 
structive or  destructive  interference  of  the  symmetric  and  antisymmetric  normal  modes 
instead  of  the  evolution  of  a single  quasi-normal  mode.  This  interferometric  nature  of  the 
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coupler  requires  two  synchronized  waveguides  with  exact  device  length  (equal  to  the  cou- 
pling length)  to  achieve  efficient  power  transfer,  thus  imposing  severe  design  and  fabrica- 
tion restrictions.  On  the  contrary,  conditions  for  the  tapered  coupler  (Eqs.  2. 1-2.3)  are  all 
inequalities  and  can  be  satisfied  by  a wide  range  of  parameters. 

The  tapered  directional  coupler  used  in  our  integration  is  somewhat  more  compli- 
cated (Figure  2.6).  The  lower  waveguide  is  an  SQW  separate  confinement  heterostructure 
(SCH)  while  the  guiding  layer  of  the  upper  waveguide  is  an  MQW  structure.  The  tapering 
of  the  coupler  is  realized  by  gradually  disordering  the  MQWs  along  the  coupler  instead  of 
changing  the  guide  thickness.  The  disordering  of  the  quantum  wells  alters  the  refractive 
index  profiles  which  are  shown  schematically  in  representative  positions  along  the  coupler 
in  Figure  2.6.  The  coupler  is  designed  in  such  a way  that  the  light  launched  into  the  SQW 
waveguide  excites  only  the  fundamental  mode  whose  power  is  concentrated  in  the  SQW 
guide  at  the  input  end  and  the  mode  evolves  adiabatically  along  the  tapered  coupler  with 
the  power  eventually  transferred  and  concentrated  in  the  MQW  guide  at  the  output  end 
[42].  The  field  profiles  of  the  normal  modes  are  also  shown  in  Figure  2.6.  At  the  input  end 
where  MQWs  are  completely  disordered,  the  fundamental  mode  field  is  concentrated  in 
the  SQW  waveguide  while  at  the  output  end  with  undisordered  MQWs,  the  field  is  con- 
centrated in  the  MQW  waveguide.  The  fabrication  of  the  tapered  coupler  using  Zn,  F,  O2 
and  N2  disordered  MQWs  have  been  investigated  in  our  group  [43].  Its  advantages,  appli- 
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Figure  2.6  The  tapered  waveguide  coupler. 
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cation  and  fabrication  in  the  laser-waveguide  integration  will  be  discussed  in  the  next  sec- 
tion. 


2.4  The  SOW  laser-MOW  waveguide  Integration 

In  the  previous  sections,  we  analyzed  the  principles  of  the  key  components  in  the 
integration  of  an  SQW  laser  and  a tapered  waveguide  coupler.  Figure  2.7  shows  the  pro- 
posed integrated  structure.  The  multilayer  structure  can  be  grown  in  a single  growth.  A 0.5 

|im  thick,  n-type  doped  GaAs  buffer  layer  is  first  grown  on  the  n'*’  GaAs  substrate,  fol- 
lowed by  an  1.5  |xm  thick,  n'  AIq  jsGao  ^fAs  lower  cladding  layer  for  the  SQW  waveguide. 
Then,  a 0.15  |im  AIq  2oGao.8oAs  lower  guiding  layer,  a 90  A GaAs  quantum  well  layer, 
and  again  a 0.15  |0.m  AIq  2oGao.8oAs  upper  guiding  layer,  all  undoped,  are  grown.  After 
growth  of  a 0.5  |xm,  undoped  AIq  35Gao  55AS  cladding  layer,  12  periods  of  65  A GaAs/190 
A AIq  4oGao  6qAs  quantum  well/barrier  layers  are  grown  to  form  the  MQW  waveguide. 
Finally,  a 0.7  |0.m  thick  AIq  35GaQ  55AS  cladding  layer  and  a 250  A GaAs  cap  layer,  both 
undoped,  are  grown.  In  this  structure,  the  SQW  and  MQW  waveguides  are  vertically 
stacked.  In  the  laser  section,  the  MQWs  are  completely  disordered  by  Zn  implantation  and 
annealing,  thus,  becoming  the  cladding  layer  of  the  SQW  laser  waveguide.  In  the  coupler 
section,  quantum  wells  are  gradually  disordered  along  the  waveguide  by  Nitrogen  implan- 
tation and  annealing  forming  a tapered  directional  coupler.  A groove  etched  between  the 
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Figure  2.7  Monolithic  integration  of  a SQW  laser  and  a tapered  waveguide 
coupler. 
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laser  section  and  the  tapered  coupler  functions  as  a partially  transmissive  mirror  for  laser 
feedback.  In  this  structure,  light  from  the  laser  is  launched  into  the  SQW  waveguide 
through  the  etched  groove  and  coupled  into  the  output  MQW  waveguide  by  another  pas- 
sive, tapered,  adiabatic  coupler. 

Such  a structure  has  unique  advantages.  As  we  pointed  out  in  Section  2.3,  the 
tapered  directional  coupler  employs  an  adiabatic  power  transfer  process  and  can  be  real- 
ized by  using  a wide  range  of  device  parameters  such  as  dimensions  and  layer  composi- 
tions for  the  two  waveguides  instead  of  a set  of  discreet  parameters  dictated  by  the 
interferometric  nature  of  the  conventional  directional  coupler.  This  allows  for  independent 
optimization  of  individual  components.  For  example,  we  used  an  SQW  SCH  waveguide 
for  the  laser  and  an  MQWs  structure  as  the  output  waveguide  for  modulation.  Moreover, 
the  vertical  configuration  of  the  two  waveguides  eliminates  the  requirements  for  compli- 
cated and  costly  material  regrowth  without  sacrificing  the  design  and  fabrication  flexibili- 
ties. 

In  the  following  chapters,  we  investigate  the  key  concepts  and  technologies 
required  for  the  SQW  laser-MQW  waveguide  integration  viz.,  Zn  disordering  of  the 
MQWs,  etched  groove  and  their  influence  on  the  threshold  of  the  SQW  laser,  and  the  asso- 
ciated fabrication  techniques. 


CHAPTER  3 

MULTIPLE  QUANTUM  WELL  DISORDERING 
3.1  Introduction 

The  laser  section  in  the  laser-coupler  integration  consists  of  a vertically  stacked 
single  quantum  well  (SQW)  waveguide  and  a multiple  quantum  well  (MQWs)  region.  The 
multiple  quantum  wells  form  an  optical  waveguide  in  the  tapered  coupler  section  and  are 
completely  disordered  in  the  laser  section  thus  becoming  a part  of  the  cladding  layer  for 
the  SQW  laser  waveguide.  The  complete  disordering  is  necessary  to  reduce  the  threshold 
current  of  the  SQW  laser.  The  characteristic  of  low  lasing  threshold  due  to  the  quantum 
effect  in  quantum  well  lasers  is  compromised  by  the  fact  that  the  quantum  well  active 
layer  that  confines  the  injected  carriers  is  so  thin  (~  100  A)  that  it  provides  little  confine- 
ment for  the  emitted  photons.  For  optical  confinement,  the  double  heterostructure  (DH) 
should  have  an  appropriate  refractive  index  difference  and  a thickness  on  the  order  of  the 
optical  wavelength,  which  is  typically  several  thousand  A for  the  GaAs/AlGaAs  material 
system.  Generally,  a separate  confinement  heterostructure  (SCH)  is  used  to  improve  the 
optical  confinement  in  quantum  well  lasers.  Figure  3.1  shows  the  A1  compositional  pro- 
files of  step  index  and  graded  index  (GRIN)  SCH  SQW  lasers.  The  thin  GaAs  quantum 
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Figure  3.1  Step  index  and  graded  index  separate  confinement  heterostructures  (SCH) 
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well  active  layer  is  sandwiched  between  two  thin,  low  A1  concentration  AlGaAs  layers  that 
serve  as  the  barrier  for  the  quantum  well  as  well  as  the  guiding  layer.  In  addition,  these  are 
surrounded  by  two  thick,  high  A1  concentration,  lower  refractive  index  AlGaAs  cladding 
layers  that  provide  the  optical  confinement. 

The  lasing  condition  can  be  expressed  as 

where 

^\E\^da 

T = -2 

oo 

J \E\^da 

— oo 

is  the  optical  confinement  factor  in  the  active  layer  a and  E is  the  optical  field  profile,  oq  is 
the  internal  loss  coefficient  in  the  laser  cavity  with  the  length  L,  and  rj  and  T2  are  the 
reflectivities  of  the  two  facet  mirrors.  The  threshold  gain  g,/,  determines  the  threshold  cur- 
rent density.  Higher  the  threshold  gain,  higher  the  threshold  current  density,  oq,  r;  and  K2 
are  mainly  determined  by  the  materials  of  the  gain  medium  and  end  facets.  This  relation 
shows  a large  optical  confinement  factor  T will  result  in  a low  threshold  current  density. 
Physically,  the  more  the  optical  field  is  concentrated  in  the  SQW  active  layer,  the  lower  the 
threshold  current  density  as  more  of  the  generated  light  experiences  the  gain  in  the  cavity. 
The  field  profile  and  the  optical  confinement  factor  is  determined  by  the  refractive  index 
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profile  of  the  structure.  Figure  3.2  shows  the  compositions  and  thicknesses  of  the  verti- 
cally stacked  MQW-SQW  structure  under  investigation.  Its  refractive  index  and  field  pro- 
files are  illustrated  in  Figure  3.3.  Note  that  the  undisordered  MQWs  region  has  a high 
refractive  index  and  the  disordering  of  the  MQWs  reduces  its  refractive  index.  Without  the 
disordering,  the  MQWs  region  itself  is  a waveguide  next  to  the  SQW  waveguide  and  the 
field  of  the  fundamental  mode  of  the  composite  SQW/MQW  structure  is  confined  in  the 
MQWs  region,  as  shown  in  Figure  3.3a,  resulting  in  a nearly  zero  confinement  of  the  light 
in  the  SQW  active  layer  of  the  laser.  Disordering  of  the  MQWs  modifies  the  refractive 
index  profile  affecting  the  field  profile  and  hence  the  optical  confinement  factor.  Only 
when  the  MQWs  are  completely  disordered  will  the  optical  field  be  well  confined  in  the 
SQW  waveguide  (Figure  3.3b).  In  this  case,  the  SQW  laser  exhibits  the  maximum  confine- 
ment factor  and  lowest  threshold  current  density. 

Disordering  the  MQWs  in  the  laser  section  is  a better  way  to  modify  the  refractive 
index  profile  than  removing  the  MQWs  by  etching.  The  latter  approach  leaves  a very  thin 
cladding  layer  for  the  SQW  laser  with  the  metal  electrode  over  it.  The  field  penetrating 
into  the  metal  as  well  as  high  doping  cladding  layer  for  ohmic  contact  would  result  in  a 
higher  loss  and  a higher  lasing  threshold.  The  disordering  of  the  quantum  well  is  per- 
formed by  introducing  impurities  into  the  quantum  well  at  high  temperature.  The  GaAs/ 
AlGaAs  heterostructures  are  in  general  very  stable.  Thermal  anneal  below  800  °C  causes 
little  intermixing.  However,  the  situation  is  completely  different  when  impurities  are  intro- 
duced in  the  structure.  Total  intermixing  i.e.,  impurity  induced  disordering  may  occur 
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Figure  3.2  The  layer  compositions  and  thicknesses  of  the  SQW/ 
MQW  structure  under  investigation. 
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mode  field  profiles  of  the  SQW/MQW  structure 
with  the  undisordered  (a)  and  disordered  (b) 
MQWs. 
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resulting  in  the  GaAs/AlGaAs  quantum  well  becoming  an  AlGaAs  alloy,  whose  A1  com- 
position is  the  average  of  the  layers. 

In  this  chapter,  we  investigate  analytically  the  dependence  of  the  lasing  threshold 
on  the  disordering  parameters.  We  will  compare  the  calculations  with  the  experimental 
results  in  Chapter  5.  The  layer  compositions  and  thicknesses  of  the  SQW/MQW  structure 
under  consideration  for  the  SQW  laser  - coupler  integration  is  illustrated  in  Figure  3.2.  We 
use  zinc  to  disorder  the  quantum  wells  as  zinc  also  functions  as  p-type  dopant  for  the  p-n 
junction  and  ohmic  contact.  Zinc  is  implanted  into  the  sample  and  then  annealed.  The 
degree  of  disordering  depends  on  the  Zn  concentration.  In  Section  3.2,  first,  we  determine 
the  Zn  profile  due  to  annealing  by  solving  Pick’s  second  law  with  a concentration  depen- 
dent diffusion  coefficient  using  finite  difference  method  (FDM).  In  Section  3.3,  with  Zn 
concentration  in  the  MQWs  region  as  a function  of  the  annealing  time  and  position,  the 
coefficient  of  Zn  induced  A1  and  Ga  interdiffusion  between  the  quantum  well  and  barrier 
are  determined  and  we  develop  a finite  difference  approach  using  the  Crank  and  Nicholson 
method  to  obtain  the  compositional  profile  of  the  disordered  MQWs.  The  equivalent 
refractive  indices  of  the  disordered  MQWs  are  then  calculated  by  solving  the  Schro- 
edinger  equation  using  the  FDM  in  Section  3.4.  In  Section  3.5,  using  the  disordering  mod- 
ified refractive  index  profiles,  the  optical  confinement  factor  in  the  laser  active  layer  is  then 
calculated  with  the  help  of  the  multilayer  stack  theory.  Finally  in  Section  3.6,  the  laser 
threshold  dependence  on  the  MQW  disordering  is  investigated. 
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3.2  Zn  Diffusion  Profiles 


The  degree  of  the  MQWs  disordering  depends  on  the  Zn  concentration.  In  this  sec- 
tion, we  model  the  Zn  profile  due  to  annealing.  Generally,  diffusion  is  the  result  of  a con- 
centration gradient  of  mobile  atoms.  The  flux  of  the  moving  atoms  is  proportional  to  the 
negative  concentration  gradient  and  expressed  by  Pick’s  law  as 

J = -DVC 

where  J is  the  flux  of  the  diffusing  atoms,  C is  the  concentration,  and  D is  the  proportional 
constant,  termed  as  the  diffusion  coefficient.  By  applying  the  continuity  condition  that  the 


rate  of  change  of  concentration  in  a given  volume,  ^ AxAyAz , must  equal  the  net  diver- 

ot 

gence  of  the  flux  of  the  atoms  flowing  in  and  out  to  the  volume,  -V-JAxAyAz,  we  get 

^ = -v.y. 

dt 

Using  Pick’s  law  and  assuming  the  diffusion  coefficient  D is  constant,  we  obtain  Pick’s 
second  law 


dt 


DV^C. 


(3.1) 


With  the  boundary  and  initial  conditions,  the  diffusion  profile  can  be  obtained  by  solving 
the  above  equation  either  analytically  or  numerically. 
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Pick’s  second  law  applies  to  most  common  diffusions  in  which  the  diffusion  coef- 
ficient remains  constant.  However,  Zn  diffusion  in  GaAs/AlGaAs  is  more  complicated  and 
does  not  have  a constant  diffusion  coefficient.  A substitutional-interstitial  diffusion  mech- 
anism has  been  suggested  to  explain  this  diffusion[44-47].  Zn,  a group  II  element,  gener- 
ally occupies  substitutionally  the  sites  of  group  III  element  Ga,  in  GaAs,  acting  as  a singly 

ionized  acceptor  s'  whose  concentration  is  S.  At  the  same  time,  a small  portion  of  Zn  exists 

interstitially,  acting  as  a singly  ionized  donor  i'*’  with  concentration  7.  The  diffusion  coeffi- 
cient of  the  substitutional  Zn  is  negligible  and  Zn  diffusion  is  due  to  the  movement  of  the 
highly  mobile  interstitial  Zn.  The  transition  of  Zn  between  the  substitutional  and  intersti- 
tial forms  are  expressed  as 

s'  + 2h  i'*’  + Vjjj 

where  h represents  holes  with  a concentration  of  p and  v^j  the  vacancies  on  the  Ga  sites 
with  concentration  V assumed  constant.  At  equilibrium,  the  law  of  mass  action  leads  to 

K^p^S  = IV. 

By  postulate,  S»I  so  we  have  p~S  for  the  material  to  be  neutral.  Thus, 

= IV. 

On  the  other  hand,  since  the  change  of  the  total  concentration  (S+T)  is  due  to  the  diffusion 


of  the  interstitials,  we  have 
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2 


dS  ^ ^ 


ar-^9r=  V 


where  D,-  is  the  diffusion  coefficient  of  the  interstitials.  Substituting  for  I using  the  last 
equation,  we  obtain 


Since  S»I,  in  effect,  the  total  concentration  C~S.  Therefore,  for  Zn  diffusion  in  GaAs/ 
AlGaAs,  we  have  a diffusion  equation  in  the  form  of  Pick’s  second  law  but  with  a concen- 
tration dependent  diffusion  coefficient,  D = Dq  C^. 

Solving  the  Zn  diffusion  equation  for  the  impurity  profile  in  the  SQW/MQW  struc- 
ture is  further  complicated  by  the  fact  that  the  medium  is  a multilayer  heterostructure. 
Since  no  closed  form  analytical  solution  is  available,  we  use  a finite  difference  method 
(FDM).  The  diffusion  equation  is  approximated  by  a group  of  finite  difference  equations 
that  relate  the  concentration  at  a point  to  those  at  the  neighboring  points  and  then  solved 
subject  to  the  boundary  and  initial  conditions.  The  approximations  are  based  on 


(3.2) 


with 


2>K^D^S 


,2 


C{xq  + Ax,0  - C(xo,0 
Ax 


(3.3) 


49 


~^C{XQ,t) 

dx 


C(xq  + Ax,t)  - 2C{xQ,t)  + C(xq  - Ax,t) 
2Ax 


(3.4) 


and 

C{x,tQ  + At)-Cix,to) 

jC(x,t,)  = . (3.5) 

In  the  FDM,  the  x-t  plane,  where  x is  the  distance  coordinate  along  the  diffusion  direction 
and  t is  the  diffusion  time,  is  divided  into  a grid  of  nodes  with  each  node  represented  by  a 
pair  of  numbers  (i,j),  as  show  in  Figure  3.4.  Thus,  Equations  (3.3),  (3.4)  and  (3.5)  can  be 
expressed  as 


. . 


Ax 


(3.6) 


2 

a c 


(Axf 


(3.7) 


and 


dt 


ij 


At 


(3.8) 


The  notation  of  i and  j is  convenient  for  computer  simulation,  especially  when  adaptive 
meshes  are  used.  For  homogeneous  materials,  substitution  of  Equations  (3.6),  (3.7)  and 
(3.8)  into  Fick’s  second  law  lead  to  algebraic  finite  difference  equations  that  can  be  solved 
for  diffusion  profiles. 
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Figure  3.4  Finite  difference  mesh  for  the  diffusion  equation. 
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For  our  SQW/MQW  multilayer  heterostructure,  the  interface  problem  should  be 
considered.  For  efficient  computing,  use  of  adaptive  meshes  is  desired.  We  derive  a set  of 
general  finite  difference  equations  that  use  adaptive  meshes  and  apply  to  both  the  points  in 
the  volume  and  at  the  interfaces.  The  Zn  diffusion  equation  (3.2)  becomes 


dC  _ d f 
dt  dx\  dx 


-3  = 


^9  C dDdC 
^^^TxTx- 


(3.9) 


We  consider  an  interface,  located  at  the  coordinate  i,  between  two  media  with  diffusion 
coefficients  and  Dy,  respectively,  as  shown  in  Figure  3.5a.  At  the  interface,  the  continu- 


ity condition  requires  that  the  flux  of  diffusing  atoms  be  continuous;  therefore. 


ac 

D 


= D 


dCb 

^dx 


= F. 


(3.10) 


For  finite  difference  expressions,  we  extend  medium  a one  step  Ax^  to  the  right  of  the 
interface  and  introduce  a ‘fictitious’  concentration  (Figure  3.5b).  Using  (3.3),  the 

left-hand  side  of  Equation  (3.10)  can  be  expressed  as 


^ 2Ax 


= F. 


(3.11) 


Also, 


(3.12) 


Eliminating  j in  (3.12)  by  using  (3.1 1)  and  substituting  into  (3.9),  we  have 
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Figure  3.5  Finite  difference  mesh  at  the  interface  i 
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. 


(3.13) 


Similarly,  using  a ‘fictitious’  concentration  Cbi.j  j (Figure  3.5c),  we  have 


. . . 


(3.14) 


Eliminating  F from  (3.12)  and  (3.13)  and  using  (3.7)  and  (3.5),  we  obtain  the  finite  differ- 
ence expression  of  the  diffusion  equation 


Therefore,  the  concentration  at  point  i and  time  j+1  can  be  calculated  from  the  concentra- 
tions of  neighboring  points  i-1,  i,  and  i+1  at  a previous  time  j.  Note  that  for  homogeneous 
materials  with  non-concentration  dependent  diffusion  coefficient  (D^=D^=constant)  and 
constant  mesh  size  (Ax^=Ax^),  the  above  equation  reduces  to  the  conventional  form  of 
finite  difference  equation  [48] 


(3.15) 


In  the  laser  section  of  the  integrated  structure,  Zn  is  first  implanted  into  the  sample 
and  then  annealed.  The  advantage  of  Zn  implantation  over  diffusion  by  Zn  vapor  or  depos- 
ited metallic  Zn  is  that  the  implant  dosage  of  Zn  can  be  easily  controlled.  When  the  anneal 
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temperature  and  time  are  fixed,  the  level  of  Zn  concentration  in  the  MQW  region  deter- 
mines the  degree  of  the  MQW  disordering  since  the  impurity  induced  interdiffusion  coef- 
ficient is  a function  of  Zn  concentration.  In  the  calculation,  we  use  an  anneal  time  of  40 
minutes  and  a temperature  of  700  °C.  We  assume  the  Zn  diffusion  coeffi- 

19  2 

cient  D{x)  = Dq{C{x)/  1x10  ) , where  Dq  is  extracted  from  references  49  and  50  for 
various  material  compositions  as  follows:  GaAs  1x10  cm  /s,  A1 35Ga  65As  IxlO  cm  / 
s,  the  MQW  5x10'*'^  cm^/s,  A1  ^Ga  gAs  8x10'^^  cm^/s.  For  the  MQWs,  an  equivalent  dif- 
fusion coefficient  is  used.The  initial  depth  of  the  implanted  Zn  is  0.35  |im.  The  layer  thick- 
nesses of  the  MQW/SQW  structure  under  consideration  are  shown  in  Figure  3.3.  Figure 
3.6  shows  the  calculated  Zn  profiles  after  annealing  with  various  initial  Zn  concentrations. 
Note  that  the  diffusion  fronts  are  all  abrupt,  which  is  the  typical  characteristic  of  Zn  diffu- 
sion in  III-V  semiconductors.  Since  the  diffusion  coefficient  increases  with  the  Zn  concen- 
tration, the  diffusion  depth  increases  with  the  initial  Zn  concentration.  The  calculated  Zn 
profile  data  as  a function  of  depth  and  time  are  saved  for  the  numerical  modeling  of  the 
MQW  disordering  in  the  next  section. 


C (cnT^) 
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Figure  3.6  Calculated  Zn  profiles  in  the  GaAs/AlGaAs  MQWs/SQW 
structure  after  annealing  with  various  implanted  Zn  concentrations  Cqo 
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3.3  Compositional  Profiles  of  Disordered  MOWs 

We  introduced  the  impurity  induced  disordering  of  quantum  wells  in  the  first  sec- 
tion (3.1)  of  this  chapter.  In  this  section,  we  discuss  the  mechanism  of  the  disordering  and 
model  the  compositional  profiles  of  the  disordered  multiple  quantum  wells  due  to  the  Zn 
diffusion,  discussed  in  the  previous  section.  The  results  will  be  used  to  investigate  the 
dependence  of  the  lasing  threshold  on  the  disordering  parameters  in  the  later  sections.  The 
impurity  induced  disordering  was  discovered  in  early  1980s  [51,  52].  Several  mechanisms 
have  been  proposed  to  explain  the  Zn  diffusion  induced  disordering  including  the  Zn 
vacancy  pair  mechanism  [51],  the  single  vacancy  ring  and  divacancy  ring  mechanism  [53], 
the  As  antisite  defect  mediated  mechanism  [54],  the  mechanism  involving  Coulombic 
interactions  between  substitutional  impurity  atoms  and  host  atoms  [55],  and  the  self-inter- 
stitial mechanism[56-59].  The  self-interstitial  mechanism  has  been  widely  accepted.  We 
discussed  in  Section  3.2  that  Zn  diffuses  in  GaAs/AlGaAs  via  the  substitutional-interstitial 
mechanism.  Zn  atoms  move  as  interstitial  donors  and  then  occupy  group  HI  vacancy  sites 
becoming  immobile  substitutional  acceptors.  The  group  HI  vacancies  are  generated  by  Zn 
interstitials  via  a ‘kick-out’  mechanism  in  which  group  III  atoms  Ga  and  A1  become  inter- 
stitial and  their  sites  become  vacancies.  Similar  to  the  diffusion  of  Zn,  Ga  and  A1  on  their 
sites  have  negligible  diffusion  coefficient  so  the  interdiffusion,  or  intermixing,  of  Ga  and 
A1  between  the  quantum  well  and  barrier  layer  is  negligible.  However,  when  Zn  is  diffused 
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into  the  crystal,  interstitial  Ga  and  A1  are  generated  via  the  kick-out  mechanism.  Like  Zn 
interstitials,  these  interstitials  also  have  large  diffusion  coefficients.  Because  of  the  con- 
centration difference  between  the  GaAs  well  (low  A1  and  high  Ga)  and  the  AlGaAs  barri- 
ers (high  A1  and  low  Ga),  the  A1  interstitials  in  the  barrier  layers  will  diffuse  into  the  well 
layers  and  the  Ga  interstitials  in  the  well  layer  diffuse  into  the  barrier  layers,  leading  to  a 
layer  disordering  (Figure  3.7).  It  is  these  Zn  diffusion  generated  A1  and  Ga  interstitials  that 
enhance  the  interdiffusion  between  the  layers. 

Pick’s  second  law  (3.1)  can  be  applied  to  the  A1  and  Ga  interdiffusion  between  the 
quantum  well  and  barrier  layers  since  it  is  a diffusion  of  interstitials.  It  is  generally 
accepted  that  the  diffusion  coefficients  for  A1  in  GaAs  and  Ga  in  AlGaAs  are  very  close 
[59].  The  interdiffusion  coefficient 


n ~ j^GaAs  _ p.  AlGaAs 

Ui  ~ U^i  ~ Uq^ 


where  D^i  and  Dq^  are  the  diffusion  coefficients  of  A1  in  the  GaAs  well  and 

Ga  in  the  AlGaAs  barriers,  respectively.  It  should  be  noted  that  the  interdiffusion  coeffi- 
cient is  Zn  concentration  dependent  according  to  the  self-interstitial  mechanism.  It  varies 
with  the  disordering  time  since  the  interdiffusion  proceeds  at  the  same  time  as  the  Zinc 
diffusion. 

Two  approaches  have  been  used  to  model  the  Zn  induced  disordering  in  the  GaAs/ 
AlGaAs  quantum  wells.  The  first  approach  uses  the  Boltzmann-Matano  technique  [60]. 
This  approach  basically  transforms  the  two- variable  diffusion  equation  complicated  by  the 
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Figure  3.7  Illustration  of  the  disordering  of  a GaAs/ 
AlGaAs  quantum  well. 
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concentration  dependent  diffusion  coefficient  into  a single-variable  equation.  This  tech- 
nique requires  that  the  boundary  conditions  can  be  described  in  term  of  the  new  variable, 
which  restricts  the  application  of  this  technique  in  many  cases.  Moreover,  the  requirement 
that  the  diffusion  coefficient  is  a function  of  concentration  only,  is  a major  problem  when 
using  this  approach  for  the  disordering  process.  As  discussed  earlier  in  this  section,  the 
interdiffusion  coefficient  varies  with  time.  Several  researchers  used  this  technique  to 
obtain  information  on  the  disordering  [59,  61].  It  has  been  pointed  out  [62]  that  the  Boltz- 
mann-Matano  analysis  does  not  apply  in  the  case  of  disordering  due  to  the  time  variance 
of  the  interdiffusion  coefficient. 

Another  approach  uses  an  error  function  approximation  [63,  64].  The  A1  composi- 
tion profile  of  the  disordered  quantum  well  is  approximated  by 


where  xq  is  the  A1  concentration  of  the  AlGaAs  barrier  before  the  disordering,  is  the 
well  thickness,  D the  interdiffusion  coefficient,  and  t the  interdiffusion  time.  Again,  this 
approach  ignores  the  time  dependence  of  the  interdififusion  coefficient.  The  assumption  for 
this  approximation  is  that  the  barrier  thickness  is  infinite,  which  is  not  the  case  for  multi- 
ple quantum  wells  and  superlattices.  Moreover,  for  slight  disordering,  the  error  function  is 
a poor  approximation. 

We  develop  a finite  difference  approach  with  a Zn-concentration  and  time  depen- 


dent interdiffusion  coefficient.  The  diffusion  equation  from  (3.1)  for  A1  is 
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^X(Z,  t)  = DiCzn,  t)-^x(z,  t) 
dz 


dt 


(3.16) 


where  A:(z,f)  is  the  A1  concentration  at  the  position  z and  time  t,  D{C2„,t)  the  interdiffusion 
coefficient  with  Zn  concentration  C2n  at  time  t.  Using  similar  arguments  for  the  Zn  diffu- 
sion, we  can  obtain  an  explicit  finite  difference  expression  of  the  diffusion  equation  identi- 
cal to  Equation  (3.15)  except  that  the  last  two  non-linear  terms  are  dropped  because  the 
interdiffusion  coefficient  D{C2„,t)  is  not  a function  of  the  A1  concentration  x. 

However,  the  explicit  expression  (3.15)  does  not  always  give  stable  and  converging 
solutions.  Care  must  be  taken  of  the  choice  of  the  mesh  sizes.  In  order  to  ensure  a stable 
solution  and  reduce  error,  we  derive  an  implicit  finite  difference  expression  using  the 
Crank  and  Nicholson  formula  [65].  When  deriving  the  explicit  expression,  values  of  x(z,t) 
at  time  j and  j+l  were  used  for  the  term  on  the  left  hand  side  of  Equation  (3.16)  while  val- 
ues at  time  j were  taken  for  the  terms  on  the  right  hand  side.  Crank  and  Nicholson  pro- 
posed that  a more  accurate  approximation  is  to  take  the  average  of  the  right  term  at  time  j 
and  at  time  j+l.  They  also  proved  that  this  approximation  always  gives  stable  solutions. 
Using  the  Crank  and  Nicholson  method  and  deriving  in  a similar  way  to  that  in  Section 
3.2,  we  obtain  an  implicit  finite  difference  expression  of  the  diffusion  equation  (3.16) 

+ + (3-17) 


where  Aj  = - 


2DAt 


(Az^  + Azb)Az^' 
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and  D is  the  interdiffusion  coefficient,  Az^  and  Az^  the  adaptive  mesh  sizes,  and  At  the  time 
interval.  For  all  the  I’s  in  the  z-t  plane.  Equation  (3. 17)  gives  an  Nt^N  matrix  equation 
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Note  that  the  boundary  conditions  have  already  been  used  in  this  equation.  If  the  concen- 
tration profile  at  time  j is  known,  then  the  concentration  at  time  y+1  can  be  calculated  from 
the  N\N  matrix  equation.  Thus,  with  the  initial  concentration  profile,  we  can  calculate  the 
profile  at  any  time.  It  should  be  noted  that  we  could  not  use  the  Crank  and  Nicholson 
method  for  the  Zn  diffusion  because  the  concentration  dependent  diffusion  coefficient 
results  in  two  non-linear  terms  in  the  equation. 
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The  Zn  concentration  profile  is  determined  in  the  previous  section  by  solving  the 
concentration  dependent  Pick’s  second  law  (3.9).  The  Zn  concentration  dependence  of  the 
A1  interdiffusion  coefficient  is  extracted  from  reference  66.  Since  the  Zn  concentration 
profile  changes  as  the  Zn  diffusion  proceeds,  we  also  obtain  a time-variant  A1  interdiffu- 
sion coefficient.  In  the  program,  we  start  to  count  the  interdiffusion  time  when  the  Zn  con- 

centration  in  the  quantum  well  under  consideration  reaches  1x10  cm  due  to  the  Zn 

diffusion  since  it  is  known  that  below  Ixio'^  cm'^,  the  interdiffusion  is  negligible.  This 
saves  the  computing  time  and  computer  memory.  The  parameters  used  in  the  calculation 
are  the  same  as  those  in  the  previous  section.  Figure  3.8  shows  the  calculated  A1  composi- 
tional profiles  of  the  MQWs  after  annealing  with  different  initial  implanted  Zn  concentra- 
tions. The  corresponding  Zn  profiles  after  the  annealing  are  shown  in  Fig.  3.6  in  the 
previous  section.  As  shown  in  Figures  3.6  and  3.8,  when  the  initial  concentration  Coq  is 

<4x10*^  cm‘^,  the  Zn  diffusion  front  is  above  the  MQWs  region  (x<0.7  |im)  and  the 
MQWs  are  not  disordered.  As  Cqq  is  increased  to  5x10*^  cm'^,  the  Zn  front  reaches  the 
MQWs.  The  Al/Ga  interdiffusion  occurs  in  the  first  quantum  well  and  other  quantum  wells 
below  the  Zn  front  remain  undisordered.  When  Cqo  is  further  increased,  the  Zn  front  goes 
deeper  into  the  MQWs  region  and  the  wells  behind  the  Zinc  front  are  gradually  disor- 
dered. At  Cqo  = 2x10^®  cm'^,  the  MQWs  are  completely  disordered  and  become  an  alloy 
of  AlGaAs  with  an  A1  composition  of  ==  0.3,  which  is  the  average  of  the  quantum  well  and 
barrier  layers. 
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Figure  3.8  Calculated  A1  compositional  profiles  in 
the  MQW  region  after  annealing  with  various  ini- 
tial implanted  Zn  concentration  Cqq 
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3.4  Refractive  Index  Profiles  of  Disordered  MOWs 

As  discussed  in  Section  3.1,  the  purpose  of  the  MQWs  disordering  in  the  laser  sec- 
tion is  to  modify  the  refractive  index  profile  of  the  MQWs/SQW  structure.  It  is  anticipated 
that  A1  and  Ga  interdiffusion  between  the  quantum  well  and  barrier  layers  induced  by  Zn 
diffusion  should  alter  the  compositional  profile  of  the  quantum  well  structure  and  cause 
changes  in  band  structure  and  energy  levels,  resulting  in  a modification  of  the  refractive 
indices  of  the  quantum  wells.  Compared  to  bulk  materials,  the  absorption  edge  of  a quan- 
tum well  structure  moves  to  higher  energy  due  to  quantized  energy  levels  above  the  bulk 
band  edge.  The  refractive  index  can  be  calculated  from  the  absorption  spectra  of  the  quan- 
tum wells.  Sophisticated  band-structure  techniques  such  as  the  pseudopotential  method 
[67],  the  tight-binding  method[68]  and  the  hybrid  band  structure  technique  [69]  have  been 
used  to  calculate  the  real  and  imaginary  parts  of  the  dielectric  constants  of  the  quantum 
wells  and  superlattices.  For  our  modeling,  we  use  an  approximation  method.The  refractive 
index  of  the  quantum  well  generally  includes  contributions  from  the  bulk  material,  the 
quantum  effect  (the  absorption  edge  shift)  and  the  exciton  effect  of  the  quantum  well.  In 
our  MQW/SQW  structure,  the  MQWs  are  designed  for  the  modulator  so  that  the  energy  of 
the  photons  emitted  from  the  SQW  laser  is  below  the  absorption  edge  of  the  MQWs.  Thus, 
in  the  absorption  spectra  of  the  MQWs,  the  overlap  between  the  exciton  peak  and  the  pho- 
ton spectra  is  small  leading  to  little  contribution  to  the  refractive  index  by  the  exciton 
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effect.  The  refractive  index  is  mainly  determined  by  the  bulk  effect  when  the  photon 
energy  is  below  the  absorption  edge  of  the  quantum  well.  This  is  the  base  for  the  average 
index  approximation  for  an  MQW  waveguide  [70].  However,  the  average  index  approxi- 
mation ignores  the  absorption  edge  shift  to  higher  energy  due  to  the  quantum  effect.  In  our 
model,  we  assume  the  refractive  index  of  the  GaAs/AlGaAs  quantum  well  at  optical 
energy  below  the  absorption  edge  to  be  equivalent  to  that  of  the  AlGaAs  alloy  that  has  the 
same  band  gap  energy.  This  approximation  has  been  used  in  the  modeling  of  other  devices 
in  our  group  leading  to  good  agreement  with  experimental  results  [42]. 

We  first  calculate  the  bandgap  of  the  disordered  quantum  well.  The  quantized 
energy  levels  in  the  conduction  and  valence  bands  can  be  calculated  by  solving  Schoe- 
dinger  equation  with  given  compositional  profiles  of  the  disordered  quantum  well.  The 
equation  can  be  written  as 

dx 

where  'F(x)  is  the  electron  or  hole  wavefunction,  V{x)  the  band  edge  potential  profile  of 
the  disordered  quantum  well  determined  by  the  compositional  profile  obtained  in  the  pre- 
vious section,  and  is  the  energy  level  to  be  solved  for.  The  boundary  conditions  for  the 
quantum  well  at  the  interface  are  determined  by  the  continuity  conditions 

Vi|.  = V2|. 


and 
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J_^l  _ J_^2 
m^dx  . m2dx 

We  use  the  finite  difference  method  to  solve  the  Schoedinger  equation.  Using  tech- 
niques similar  to  those  presented  in  Sections  3.2  and  3.3,  we  derive  a general  finite  differ- 
ence expression  for  Schroedinger  equation  accommodating  the  interface  continuity 
conditions  and  the  adaptive  meshes.  The  finite  difference  expression  for  the  Schoedinger 
equation  can  be  written  as 


where 


.4  - ^ 


(3.18) 


Ax 


2 

^ . ( — ]—  + — (3.19) 

,1 -H  Ax, .2^,1  Ax., 


r 

m,.2Ax,-2(Ax.j  -I-  Ax, 2) 


(3.20) 


and  is  the  potential  value  at  point  i,  E is  the  energy,  and  Ax,-y  and  Ax, -2  the  adaptive  mesh 
sizes  at  point  /.  For  all  the  i’s,  the  finite  difference  equation  becomes  an  NxN  matrix  equa- 
tion 
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The  boundary  condition  that  \|/o  and  \|/n+j  are  zero,  has  been  used.  In  the  matrix  equation, 
for  \|/i,  \|/2v,  tj/n-i  Vn  have  non-trivial  solution,  the  determinant  of  the  NxN  matrix 
must  be  zero. 


5,  C, 

A2  B2  C2 


= 0. 


^n-l  ^n-]  ^n-1 


(3.21) 


Note  that  A,-,  5,,  and  Ci,  i=l,  2,...,  N in  Equations  (3.18),  (3.19),  and  (3.20)  are  functions  of 
unknown  energy  E.  Solving  Equation  (3.21)  for  E gives  a set  of  discreet  energy  values  E„ 
which  are  the  allowed  energy  levels  in  the  quantum  well. 

Once  the  energy  levels  in  the  conduction  and  valence  bands  and  hence  the  bandgap 
of  the  disordered  quantum  well  are  obtained,  we  can  calculate  the  refractive  index  using 
Afromowitz’s  refractive  index  model  for  the  Al^Caj.^As  material  [71],  according  to  our 
approximation  approach.  In  the  calculation,  we  assume  the  parabolic  band  structure  and  a 
ratio  of  0.62:0.38  is  used  for  barriers  in  the  conduction  band  and  the  valence  band.  The 
potential  profile  V{x)  or  V}  of  the  disordered  quantum  well  is  obtained  from  the  composi- 
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Figure  3.9  Calculated  refractive  index  profiles  of  the  MQW/ 
SQW  structure  with  the  MQWs  disordered  under  various 
initial  Zn  concentration  Cqq 
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tional  profile  calculated  in  the  previous  section.  Figure  3.9  shows  the  calculated  refractive 
index  profiles  of  the  MQW/SQW  structure  with  the  MQWs  disordered  under  various  ini- 
tial implanted  Zn  concentrations.  The  corresponding  Zn  concentration  profiles  and  the  A1 
compositional  profiles  of  the  disordered  quantum  wells  are  shown  in  Figures  3.6  and  3.8. 
As  shown  in  Figure  3.9,  the  refractive  index  of  the  MQWs  are  gradually  modified  along 
the  depth  direction  as  the  initial  implanted  Zn  concentration  is  increased  and  the  Zn  front 

>j(\  -5 

goes  deeper.  When  the  initial  concentration  is  increased  to  ~ 2x10  cm  , the  MQWs  are 
completely  disordered  and  its  refractive  index  reduces  to  that  of  an  Al^Caj.^As  alloy  with 
x=0.30.  This  would  dramatically  improve  the  optical  confinement  of  the  lasing  light  in  the 
SQW  active  layer. 


3.5  The  Optical  Confinement  Factor 

As  we  discussed  in  Section  3.1,  the  threshold  current  density  of  a laser  is  inversely 
proportional  to  the  optical  confinement  factor  which  is  defined  as 

j\E\^da 

r = -2 

oo 

J \Ef-da 

— oo 

where  E is  the  field  profile  of  the  normal  mode  in  the  SQW/MQW  structure  and  a repre- 
sents the  SQW  active  layer.  When  the  refractive  index  profile  of  the  MQWs  is  modified  by 


70 


the  Zn  induced  disordering,  the  field  profile  is  altered  accordingly  leading  to  a change  in 
the  optical  confinement  factor.  We  calculate  the  confinement  factors  with  disordered 
MQWs  modeled  in  the  previous  sections. 

Although  the  MQWs  and  the  SQW  SCH  in  the  MQW/SQW  structure  are  two  sin- 
gle mode  waveguides,  the  whole  structure  can  be  considered  a multilayer  waveguide,  as 
we  discussed  in  Sections  2.1  and  2.3.  The  normal  mode  field  profile  of  the  multilayer 
waveguide,  representing  the  optical  field  distribution  in  the  MQW/SQW  structure,  can  be 
calculated  using  the  multilayer  stack  theory.  Figure  3.10  shows  a ^-layer  waveguide  with 
the  coordinates  in  the  transverse  direction  of  the  waveguide.  From  the  wave  equation,  the 
field  profile  in  each  layer  can  be  expressed  as 


N <n- 
N>ri; 


i = 1,2,  ...,q 


(3.22) 


where 


u-  = K-{x-x^), 


N = ^ , the  mode  index, 
ko 


and  P is  the  propagation  constant,  and  n,  the  refractive  index  in  layer  i.  ai  and  are  the 
coefficients  determined  by  the  boundary  conditions.  At  the  interface  between  layer  i and 
i+l,  the  field  must  satisfy  the  boundary  conditions 
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Figure  3.10  A ^-layer  waveguide. 
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and 

These  two  boundary  conditions  lead  to  a relationship  between  the  field  amplitudes  in  lay- 
ers i and i+l 


MAN) 


*1  + 1 


b; 


i+l 


(3.23) 


where  is  the  characteristic  matrix  of  the  layer  whose  elements  are  functions  of  the 
mode  index  N.  Physically,  for  a guided  mode,  the  field  must  exponentially  decay  in  the 
outer-most  layers  (i=l  and  q).  Since  in  these  two  layers  A^>n,-,  the  fields  are  expressed  in 
terms  of  cosh  and  sink  functions  in  (3.22).  To  represent  the  exponentially  decaying  fields 
in  these  two  layers,  the  coefficients  in  (3.22)  should  satisfy  aj  = bj  and  =-  b^,  so  that 
under  the  coordinate  system  shown  in  Figure  3.10, 


and 


£j(x)  = UjCoshM,  + sinhMj 


fljg 


1 


= a^coshw^  + ^j^sinhM^  = a^e 


On  the  hand,  if  we  use  Equation  (3.22)  successively  from  layer  7 to  layer  q,  we  can  get 
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a. 


= Mi(AT)M2(iV)...M  ,(iV) 


(3.24) 


If  we  use  the  decaying  conditions,  i.e.,  aj=bi  and  a^=-bg,  and  assign  an  arbitrary  value  to 
a^,  say,  a^=l,  then  Equation  (3.24)  becomes 


a. 


1 

-1 


(3.25) 


The  two  equations  in  (3.25)  can  then  be  solved  for  the  two  unknowns  N and  a j.  The  values 
of  N are  the  effective  indices  of  the  normal  modes.  Once  N and  O;  are  obtained,  the  field 
profiles  of  the  normal  modes  can  be  calculated  using  Equations  (3.23)  and  (3.22).  Note 
that  we  assigned  an  arbitrary  value  to  so  the  field  profile  is  not  normalized.  But  this 
does  not  influence  the  results  of  the  confinement  calculation. 

Figure  3.1 1 shows  the  calculated  optical  field  profiles  in  the  laser  section  with  dis- 
ordered MQWs.  The  refractive  index  profiles  of  the  disordered  MQWs  are  obtained  from 

the  previous  sections.  At  Cqq  = 4x10*^  cm'^,  the  MQWs  are  not  disordered  and  have 
higher  refractive  indices  than  those  of  the  surrounding  layers,  forming  a guiding  layer  with 
the  field  of  the  fundamental  mode  well  confined  in  it.  This  results  in  an  optical  confine- 
ment factor  of  = 0 in  the  SQW  active  layer.  When  Cqo  increases,  the  MQWs  are  gradually 
disordered  and  the  refractive  index  profile  changes  accordingly,  resulting  in  an  alteration 
in  the  field  profile  and  the  optical  confinement  factor  in  the  SQW  layer.  At  Cqq  = 9.0x10^^ 
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Figure  3.1 1 Calculated  optical  field  profiles  in  the  SQW/MQWs 
structure  with  the  MQWs  disordered  under  various  initial  Zn 
concentration  Cqo- 
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cm’^,  the  field  in  the  SQW  waveguide  starts  to  increase  but  most  of  the  power  remains  in 
the  MQW  waveguide  leading  to  a small  optical  confinement  factor  in  the  SQW  guide.  At 

Cqo  = 9.8x10^^  and  9.9x10^^  cm'^,  the  field  is  roughly  even-split  in  the  MQW  and  SQW 

waveguides.  At  Cqq  = 1.1x10^®  cm'^,  most  of  the  field  switches  over  to  the  SQW 

waveguide.  The  quick  change-over  of  the  field  distribution  around  Coo=9.8xlO  cm"  due 
to  the  MQW  disordering  leads  to  a rapid  increase  in  the  optical  confinement  factor  in  the 
SQW  active  layer.  At  Cqo  = 1.1x10^®  cm"^,  although  the  MQWs  are  not  completely  disor- 
dered yet,  the  field  is  well  confined  in  the  SQW  waveguide  and  the  confinement  factor  is 
close  to  the  maximum.  At  Cqq  = 2x10^®  cm"^,  the  MQWs  are  completely  disordered  and 
become  an  alloy  AlGaAs  with  an  A1  concentration  of  « 0.3  whose  refractive  index  is 
slightly  higher  than  those  of  its  cladding  layers.  The  field  profile  does  not  change  much 

from  that  for  the  case  of  Cqq  = 1.1x10^®  cm"^.  Figure  3.12  shows  the  variation  of  the  opti- 
cal confinement  factor  with  Coq.  As  we  analyzed  earlier,  the  optical  confinement  factor 

increases  sharply  at  around  Cqo  = 9.8x10*^  cm"^  and  remains  nearly  constant  for  Cqo  equal 
to  and  greater  than  1.1x10^®  cm'^. 


Optical  confinement  factor 
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Figure  3.12  Calculated  optical  confinement  factor  vs.  Cqo. 
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3.6  The  Lasing  Threshold  Current  Density 

We  have  calculated  the  optical  confinement  factor  for  the  light  in  the  SQW  active 
layer  of  the  MQW/SQW  structure  with  different  degree  of  the  MQW  disordering.  In  this 
section,  we  investigate  the  variation  of  the  lasing  threshold  current  density  with  the  disor- 
dering parameters.  The  lasing  condition  can  be  expressed  as 

+ (3-26) 

where  g,;,  is  the  threshold  gain,  F is  the  optical  confinement  factor,  Oj  is  the  internal  loss  in 
the  laser  cavity,  L is  the  cavity  length,  and  rj  and  K2  are  the  reflectivities  of  the  two  reflec- 
tion facets.  The  terms  on  the  right-hand  side  of  Equation  (3.26)  are  basically  determined 
by  the  material,  its  quality,  and  the  device  dimensions.  With  the  conflnement  factor  calcu- 
lated from  the  previous  section,  we  can  determine  the  threshold  gain  g,;,  using  Equation 

(3.26).  Through  the  relation  between  the  gain  and  the  injected  carrier  density,  we  can 
obtain  the  threshold  current  density. 

We  first  outline  the  gain  calculation.  As  we  discussed  in  Section  2.2,  in  a quantum 
well,  the  energy  band  splits  into  a series  of  subbands.  The  transition  rate  of  an  electron 
between  a conduction  band  state  and  a valence  band  state  is  given  by  Fermi’s  golden  rule 

[72]  and  it  is  proportional  to  the  transition  matrix  element  \Mj]^  which  is  determined  by 
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the  band  structure  and  the  transition  energy  Egf,  = hoi  For  a given  transition  energy  hoa,  the 
density  of  the  conduction  and  valence  state  pairs  is  given  by  a reduced  density  of  states 
Pred  which  is  determined  by  the  density  of  states  in  the  conduction  band  and  valence  band. 
Thus,  the  total  transition  rate  W is  determined  by  the  product  of  the  transition  matrix  ele- 
ment, the  reduced  density  of  states  and  the  Fermi-Dirac  distributions  and  /^,  in  the  con- 
duction and  valence  bands,  the  latter  giving  the  probability  of  a state  being  occupied  by  an 
electron.  A downward  transition  from  the  conduction  band  to  the  valence  band  generates  a 
photon  while  an  upward  transition  from  the  valence  band  to  the  conduction  band  absorbs 
one.  The  optical  gain  g is  defined  as  the  fractional  increase  in  photons  per  unit  length, 
which  is  determined  by  the  difference  of  the  downward  and  upward  transition  rates  and 
can  be  expressed  as 

g(/ico)  = (3.26) 

For  detailed  derivations  and  expressions,  the  readers  are  referred  to  reference  35. 

To  determine  the  Fermi  levels  for  the  Fermi-Dirac  distributions,  we  use  the  rela- 
tions 

oo  oo 

N = 2\f>J,dE  = J f^dE  (3.27) 

Ec  " E, 


and 
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P = 2 j p,(  1 - f,)dE  = 2p,X  J ( ^ (3-28) 

" -oo 

where  N and  P are  injected  electron  and  hole  densities,  respectively,  and  and  are 
Fermi  distributions  with  unknown  quasi-Fermi  level  % and  Efv  Efc  and  can  be  numer- 
ically calculated  from  Equations  (3.27)  and  (3.28)  for  given  N and  P.  Note  that  we  used  the 
fact  that  for  quantum  wells,  the  density  of  states  is  constant  and  the  energy  bands  are 
divided  into  subbands. 

From  Equation  (3.26),  we  obtain  the  gain  spectra  for  given  densities  of  injected 
carriers.  Because  of  the  spectral  broadening  associated  with  each  transition  of  energy  hco, 
a convolution  of  g with  a Lorentzian  lineshape  function  is  performed.  The  calculated  max- 
imum gain  in  the  gain  spectrum  versus  the  injected  carrier  density  N is  plotted  in  Fig- 
ure 3.13.  The  calculation  is  based  on  the  90  A GaAs/Alo  2Gao.8^s  single  quantum  well 
laser  used  in  the  integration.  Under  the  steady  state  condition,  the  injected  current  density 
J and  the  carrier  density  N are  related  by 

J = (3.29) 

where  is  the  well  thickness  and  the  carrier  lifetime  %=\I{BN).  Using  (3.29),  the  g,^  ~ N 
curve  (Figure  3.13)  is  converted  to  the  gmax~J  shown  in  Figure  3.14. 

The  threshold  gain  is  determined  by  (3.25)  with  the  optical  confinement  factor  cal- 
culated from  the  previous  sections.The  corresponding  threshold  current  density  then  can 
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Figure  3.13  Calculated  peak  gain  g^j^x  io  the  gain 
spectrum  versus  the  injected  carrier  density  N. 
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J (A/cnr^) 


Figure  3.14  Calculated  as  a function  of  injected  current 
density  J 


(mA) 
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Figure  3.15  Calculated  Ijj,  vs.  Coq.  (50p,m  x 500  ^im 
stripe) 
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be  obtained  from  the  gain  calculation  outlined  above.  Figure  3.15  shows  the  calculated 
threshold  current  density  as  a function  of  the  initial  Zn  concentration.  The  threshold  cur- 
rent density  decreases  rapidly  as  the  initial  concentration  of  implanted  Zn  increases  and 
reaches  the  minimum  when  the  MQWs  are  completely  disordered.  The  sharp  decrease  in 
the  threshold  current  density  corresponds  to  the  sharp  increase  in  the  optical  confinement 
factor,  which  is  the  result  of  the  refractive  index  decrease  of  the  cladding  layer  by  the 
MQW  disordering. 

In  this  chapter,  we  modeled  the  Zn  diffusion  induced  disordering  in  the  MQW/ 
SQW  structure  and  investigated  the  variation  of  the  threshold  current  density  of  the  SQW 
laser  with  the  MQW  disordering.  We  predicted  a rapid  reduction  in  the  lasing  threshold 
when  the  MQW  cladding  layer  of  the  SQW  laser  is  totally  disordered.  The  results  of  the 
numerical  simulation  will  be  compared  with  the  experimental  results  in  Chapter  5. 


CHAPTER  4 

GROOVE  BETWEEN  THE  LASER  AND  THE  TAPER 


The  Zn  disordered  MQW  region  discussed  in  Chapter  3 acts  as  the  p contact  clad- 
ding layer  for  the  laser  section.  However  we  still  need  to  provide  a mechanism  for  forming 
the  laser  cavity  as  well  as  a means  of  coupling  the  laser  output  to  the  tapered  waveguide 
section.  The  groove,  etched  between  the  laser  section  and  the  tapered  coupler  (Figure  2.7), 
plays  a key  role  in  this  connection.  It  functions  as  a partially  transmissive  mirror.  The 
reflection  from  the  groove  provides  the  necessary  feedback  for  the  lasing  action  and  the 
transmission  couples  light  from  the  laser  to  the  tapered  coupler.  In  this  chapter,  we  investi- 
gate the  reflectance  from  the  groove  and  its  influence  on  the  lasing  threshold  and  power 
output. 


4.1  The  Characteristic  Matrix 

In  this  section,  we  first  discuss  the  reflection  and  transmission  from  a three  layer 
structure  using  a matrix  method  [73]  and  then  apply  the  matrix  method  to  the  reflection 
problem  for  the  groove  in  the  next  section. 

The  three  layer  structure,  with  a layer  of  thickness  d and  refractive  index  rij,  sand- 
wiched between  two  infinite  thick  layers  of  refractive  indices  hq  and  ri2,  is  shown  in  Figure 
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Figure  4.1  Reflection  and  transmission  from  a 
three  layer  structure. 
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4.1.  The  electric  field  incident  at  bulk  interfaces  a and  b between  the  three  layers  is 
reflected  from  and  transmitted  through  the  interfaces.  The  field  in  each  region  can  be 
expressed  as  two  waves:  forward  traveling  and  backward  traveling,  as  shown  in  Figure  4.1. 
If  light  is  launched  into  the  three-layer  structure  from  region  0,  the  forward  traveling  wave 
represents  the  transmitted  light  while  the  backward  traveling  wave  represents  the  reflected 
light.  The  notation  used  to  designate  the  waves  is  as  follows.  The  number  in  the  subscript 
represents  the  region  where  the  wave  travels,  the  letter  represents  the  interface  on  which 
the  wave  is  incident,  and  the  + and  - signs  represent  forward  and  backward  traveling 
waves,  respectively.  The  total  electric  fields  at  the  interface  a in  region  0 and  1 thus  can  be 
expressed  as 

^Oa  = ^Oa 

and 


For  normal  incident,  the  field  is  continuous  at  the  interface,  thus 


(4.1) 


For  magnetic  fields  at  the  interface,  similar  relations  hold  and  we  have 


yo(^Oa  = yi(^la 


(4.2) 


where  yg  and  yj  are  the  wave  admittance  in  regions  1 and  2,  respectively.  From  Equations 
(4.1)  and  (4.2),  we  can  express  the  forward  and  backward  traveling  waves  in  the  region  0 
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in  terms  of  those  in  region  1 in  a matrix  form. 
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where 


(4.3) 
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and 


“ «o  + «i’ 

When  the  forward  and  backward  waves  travel  between  the  two  interfaces  a and  h,  the 

accumulated  phases  and  the  fields  at  the  two  interfaces  are  related  by 

+ + -fSi 

- ^\b  ^ 

and 

where  6j  = kQti^d  is  the  phase  change.  They  can  be  expressed  in  a matrix  form 


+ 

II 

e 0 

^ib 

^ 7'5i 



1 

JP 
1 

0 e [ 

(4.4) 


Using  Equations  (4.3)  and  (4.4),  we  can  relate  the  fields  in  the  regions  0 and  2, 
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Therefore,  the  reflectance  R and  transmittance  T for  the  three  layer  structure  can  be  calcu- 
lated from  this  equation.  We  call  the  matrices  in  Equations  (4.3)  and  (4.4)  the  characteris- 
tic matrices.  By  successively  using  the  characteristic  matrices,  we  can  obtain  the 
reflectance  and  transmittance  of  a multilayer  structure. 

4.2  Reflectance  from  the  Groove 

In  semiconductor  lasers,  usually  two  cleaved  facets  are  used  as  reflectors  to  form  a 
laser  cavity.  The  reflectance  from  the  facets  is  determined  by  the  crystal  material.  For 
GaAs/AlGaAs  lasers,  the  reflectance  of  the  cleaved  facets  is  about  0.3. 

In  our  integrated  structure,  a groove  etched  between  the  laser  section  and  the 
tapered  coupler  functions  as  a reflector  while  the  other  reflector  is  formed  by  a cleaved 
facet.  The  groove,  together  with  the  laser  and  the  tapered  coupler,  forms  a multilayer 
structure  shown  in  Figure  4.2.  We  can  calculate  the  reflectance  of  the  groove  using  the 
matrix  method  discussed  in  the  previous  section.  Note  that  the  refractive  index  in  the 
tapered  waveguide  coupler  is  not  constant  along  the  taper  because  of  the  tapering.  We  cal- 
culated the  mode  index  profile  of  the  fundamental  mode  along  the  taper  using  the  Multi- 
layer stack  theory  discussed  in  Section  3.5.  The  index  profile  along  the  taper  is  shown  in 
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Figure  4.2  The  multilayer  model  for  the  reflectance  calcu- 
lation. The  taper  section  is  approximated  by  a series  of 
small  steps. 
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Figure  4.3.  The  tapered  waveguide  coupler  is  designed  so  that  the  evolution  of  the  field 
propagating  along  the  taper  remains  adiabatic  and  the  light  launch  into  the  taper  as  the 
fundamental  mode  remains  in  the  fundamental  mode  as  it  propagates  along  the  tapered 
coupler. 

The  reflectance  due  to  the  groove  is  defined  as 


el 

where  £l'*’  is  the  forward  traveling  wave  incident  on  the  groove  and  8l'  is  the  backward 
traveling  wave  reflected  from  the  groove,  as  shown  in  Figure  4.2.  In  the  calculation  out- 
lined in  the  last  section,  using  the  multilayer  matrix  method,  the  mode  index  of  the  funda- 
mental mode  in  the  MQW/SQW  structure  with  MQWs  completely  disordered  is  used  as 
the  refractive  index  of  the  laser  section  (layer  0 in  Figure  4.2).  The  tapered  waveguide  cou- 
pler is  divided  into  many  small  sections  along  the  taper  direction  so  that  the  mode  index 
profile  in  Figure  4.3  is  approximated  by  a series  of  small  steps.  The  reflectance  of  the 
groove  is  calculated  by  repetitively  using  the  characteristic  matrices  defined  in  the  last  sec- 
tion. The  calculated  reflectance  versus  the  groove  width  is  shown  in  Figure  4.4.  The  reflec- 
tance oscillates  with  respect  to  the  groove  width  since  the  groove  forms  a Fabry-Perot 
cavity.  It  should  be  noted  in  the  figure  that  the  reflectance  of  the  groove  can  be  signifi- 
cantly higher  than  that  of  a cleaved  facet  (-0.3).  This  should  reduce  the  lasing  threshold. 
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Figure  4.3  The  mode  indices  along  the  tapered  waveguide 
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Figure  4.4  Calculated  reflectance  of  the  groove  versus  the 
groove  width. 
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4.3  The  Lasing  Threshold  Variation 


The  relationship  between  the  lasing  threshold  and  the  reflectance  of  the  reflectors 
can  be  expressed  as 


where  g,;,  is  the  threshold  gain,  F the  optical  confinement  factor,  cq  the  internal  loss,  L the 
cavity  length,  and  rj  and  r2  the  reflectance  of  the  two  facets.  This  relation  shows  higher 
reflectance  results  in  lower  lasing  threshold  as  more  light  is  reflected  so  that  the  cavity 
experiences  less  external  loss. 

To  analyze  the  variation  of  the  threshold  current  density  with  the  groove  width,  we 
use  a simple  approach,  instead  of  the  rigorous  analysis  in  Chapter  3.  Generally,  the  thresh- 
old current  density  J,/,  is  proportional  to  the  threshold  gain  g,/,,  7,;,  = From  Equa- 
tion (4.5),  we  then  have 


The  threshold  current  densities  of  the  two  lasers  from  the  same  wafer,  one  with  an  etched 
groove  and  the  other  with  two  cleaved  facets  (Jth)  can  be  related  as 


(4.5) 
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(4.6) 


where  Tq  and  Tg  are  reflectivities  of  the  cleaved  facet  and  the  groove,  respectively.  With  an 

assumed  J,h  we  can  calculate  the  variation  of  7,;,^  with  the  groove  width  d from  Equation 
(4.6)  since  the  r„  vs.d  relation  is  known  from  the  calculation  in  the  previous  section.  We 

o 

assume  =4(X)  mA.  The  calculated  threshold  current  of  the  laser  as  a function  of  the 
groove  width  is  shown  in  Figure  4.5. 


In  this  section,  we  investigate  how  the  groove  width  affects  the  ratio  between  the 
power  output  from  the  taper  side  (Pj)  and  from  the  laser  side  (Pl)  (Figure  4.6).  On  the 
laser  side,  the  reflection  is  provided  by  a cleaved  facet  with  a reflectance  Rl  and  transmit- 
tance Tl-  On  the  taper  side,  the  feedback  for  the  laser  is  provided  by  the  composite 
groove-taper  structure  (Figure  4.2)  which  is  represented  by  equivalent  reflectance  R-j-  and 
transmittance  Tj  (Figure  4.6)  calculated  in  Section  4.2.  The  forward  and  backward  travel- 
ing waves  in  each  region  are  also  shown  in  Figure  4.6.  The  power  output  from  each  side 
can  thus  be  expressed  as 


4.4  Power  Output  Variation 
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Figure  4.5  Calculated  threshold  current  versus  the 
groove  width. 
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Figure  4.6  Fields  and  output  power  in  the  laser  cavity. 


Figure  4.7  Calculated  ratio  of  power  output  from  the  taper  side 
Py  and  from  the  laser  side  Pl  vs.  the  groove  width. 
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and 


In  the  laser  cavity  with  net  gain,  g-a,  the  fields  at  the  two  facets  are  related  by  a 
factor  or  depending  on  the  direction  of  wave  propagation.  With  the  lasing 

condition 

the  ratio  between  the  power  output  from  the  two  sides  is  found  to  be 

^ _ ^tJ^l 

Since  R-j-  and  Tj  are  functions  of  the  groove  width,  the  power  ratio  varies  with  the  groove 
width.  Figure  4.7  shows  the  calculated  power  ratio  as  a function  of  the  groove  width.  The 
ratio  varies  between  0.7  and  1.7. 

In  this  chapter,  we  have  investigated  the  dependence  of  the  reflectance  from  the 
groove  and  the  corresponding  lasing  threshold  current  on  the  groove  width.  As  expected, 
both  the  reflectance  and  the  threshold  current  oscillate  with  the  groove  width.  Potentially, 
their  performance  can  be  significantly  improved  over  their  counterparts  viz.,  the  cleaved 


facet  lasers. 


CHAPTER  5 

FABRICATION  TECHNIQUES  AND  EXPERIMENTAL  RESULTS 


In  this  chapter,  following  the  discussion  on  material  growth,  we  investigate  the  key 
fabrication  techniques  and  the  associated  experimental  results  including  the  MQW  disor- 
dering by  Zn  implantation  and  annealing,  the  groove  etching  by  reactive  ion  etching,  and 
the  fabrication  of  the  tapered  waveguide  coupler.  We  conclude  the  chapter  with  the  experi- 
mental results  on  the  integration  of  the  SQW  laser  and  the  tapered  waveguide  coupler. 

5.1  Material  Growth 


All  epitaxial  growths  for  our  samples  are  performed  by  metalorganic  chemical 
vapor  deposition  (MOCVD),  also  known  as  metalorganic  vapor  phase  epitaxy  (MOVPE). 
Tri-Methyl  Aluminum  (TMA),  Tri-Methyl  Gallium  (TMG),  and  Arsine  (ASH3)  are  used  as 
the  precursors.  The  metalorganic  sources,  TMA  and  TMG,  are  liquid  at  room  temperature. 
Hydrogen  carrier  gas  is  bubbled  to  vaporize  them  and  the  vapor  is  then  led  to  the  reaction 
tube.  ASH3  is  gas  at  room  temperature  and  its  flow  to  the  reactor  is  controlled  by  valves. 
The  gases  are  pyrolized  on  the  hot  surface  of  the  substrate,  generally  at  500  - 800  °C 
heated  by  RE  induction.  The  pyrolysis  reaction  is  non-reciprocal. The  crystal  grows  contin- 
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uously  only  on  the  hot  susceptor.  The  growth  at  other  places  is  suppressed  by  the  water 
cooled  tubes.  Since  the  gas  flow  rates  in  the  reactor  can  be  easily  controlled,  the  composi- 
tions of  the  alloys  and  doping  levels  can  be  changed  quickly  leading  to  narrow  transient 
regions  and  abrupt  hetero-interfaces.  MOCVD,  as  well  as  MBE,  makes  the  growth  of 
quantum  well  and  superlattices  structure  possible. 

The  substrate  is  cleaned  and  etched  before  growth.  The  procedure  is  as  follows, 

1)  Dip  in  boiling  Trichloroethane  (TCE)  for  5 minutes, 

2)  dip  in  boiling  Aceton  for  5 minutes, 

3)  dip  in  warm  methanol  for  2 minutes, 

4)  rinse  in  DI  water, 

5)  etch  in  the  NH40H:H202:H20  solution  with  ratio  of  20:7;  100  for  5 minutes 

to  remove  possible  damage  on  the  wafer  surface  caused  by  lapping  and  polishing 
of  the  wafer, 

6)  rinse  in  DI  water  and 

7)  blow  dry  by  N2. 

We  use  an  atmospheric  pressure  horizontal  reactor.  The  growth  is  at  a temperature 
of  725  °C.  The  growth  rate  used  is  300  A/minute  with  a V/III  mole  ratio  of  40.  n-type 
doped  GaAs  wafers  of  (100)  orientation  are  used  as  substrates.  Figure  5.1  shows  the  P-I 
curve  of  a SQW  laser  (without  the  MQWs)  we  grew.  In  addition  to  the  samples  we  grew 
here  at  UF,  we  also  used  the  samples  from  outside  vendors. 
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Figure  5.1  P-I  curve  of  a regular  SQW  laser  ( 1 00  |.im  x 500  |im). 
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s 9 MOW  Disordering  bv  Zn  Implantation  and  Annealum 

In  Chapter  3,  we  modeled  the  Zn  induced  disordering  of  the  MQWs.  Zn  is  first 
introduced  into  the  sample  by  ion  implantation  and  the  subsequent  annealing  disorders  the 
MQWs.  Ion  implantation  is  the  introduction  of  energetic  charged  particles  into  targets 
with  sufficient  energy  to  penetrate  beyond  the  surface  region.  A beam  of  dopant  ions  with 
fixed  energy  is  rastered  across  the  sample  surface.  The  penetration  depth  is  detenmned  by 
the  energy  of  the  incident  ions.  They  enter  the  crystal  lattice,  collide  with  the  atoms  in  the 
lattice,  and  gradually  lose  their  energy  coming  to  rest  at  some  depth  in  the  crystal.  The 
dose  of  the  implanted  ions  is  determined  by  the  ion  flux  and  the  implantation  time.  The 
flux  of  the  charged  ions  forms  a current  that  can  be  directly  measured.  Thus,  the  dose  of 
the  implantation  can  be  easily  controlled.  Also,  high  concentration  doping  can  be  achieved 
with  implantation.  For  these  reasons,  we  choose  implantation  over  other  methods  such  as 
diffusion  from  vapor  or  deposited  Zn. 

The  Zn  implantations  in  samples  with  the  MQW/SQW  structure  (Figure  3.3)  have 
been  performed  with  various  dose  levels.  The  implantation  energy  has  been  kept  at  165 
keV  which  is  the  highest  energy  attainable  from  the  implanter  we  used.  After  the  Zn 
implantation,  the  sample  is  capped  by  500  A thick  Si02  deposited  by  e-beam  evaporation. 
The  sample  is  then  placed  top-side  down  on  a silicon  wafer  and  annealed  in  a furnace  at 
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700  °C  for  40  minutes.  The  Si02  cap  is  used  to  prevent  the  loss  of  asenic  from  the  sample 
surface  during  the  annealing. 

The  SQW  laser  with  the  MQW  cladding  layer  is  then  fabricated  from  the  Zn 
implanted  and  annealed  sample.  The  procedure  is  as  follows. 

1)  Remove  the  Si02  cap  layer  by  Buffered  Oxide  Etchant  (BOE)  6:1,  etching  for 

10  seconds,  then  rinse  in  D1  water  and  blow  dry  with  N2. 

2)  Spin  photoresist  S 1400-33  at  5000  RPM  for  40  sec.  Bake  at  90  °C  for  30  min. 
The  thickness  of  the  photoresist  is  about  1.5  p,m. 

3)  Expose  on  the  mask  aligner  for  7 sec.  The  mask  pattern  is  50  p.m  wide  stripes 
for  laser  electrodes. 

4)  Develop  in  MF-319  for  40  sec.  Rinse  in  DI  water  and  Blow  dry  with  N2. 

5)  E-beam  evaporation  of  Ti  (200  A)  and  Au  (1000  A). 

6)  Lift-off  in  aceton  to  get  50  |im  wide  Ti/Au  electrodes. 

7)  Lap  the  backside  until  sample  is  -100  |xm  thick. 

8)  Clean  the  sample  in  warm  TCE,  aceton  and  methanol.  Rinse  in  DI  water.  Blow 
dry. 

9)  E-beam  evaporate  of  AuGe  for  1000  A on  the  backside. 

10)  Cleave  the  sample  into  500  |xm  laser  bars  for  test. 

We  used  various  dosage  of  Zn  implantation.  At  low  doses  (below  10^^  cm'^)  no 
lasing  was  observed.  When  the  dose  increased  to  the  order  of  10*^,  the  device  lased  with 


Power  P (a.  u.) 
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Figure  5.2  P-I  curve  of  a MQW/SQW  laser  with  a implant 
dosage  of  7xl0'^  cm'^. 
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Figure  5.3  vs.  Coq  for  ttie  .SQW  laser  (50  pm  x 500  |im) 
with  the  disordered  MQW  clading.The  dots:  experimental. 
The  solid  line:  calculated. 
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high  threshold  current  indicating  partial  disordering  of  the  MQWs.  When  the  dose  was 
increased  from  3x10*^  cm'^  to  7x10*^  cm’^,  which  correspond  to  initial  concentration  Coo 

of  implanted  Zn  of  8x10*^  to  2x10^®  cm’^,  the  threshold  current  decreased  rapidly  from 
-1000  mA  to  -400  mA.  Figure  5.2  shows  the  P-I  curve  of  the  laser  with  a Zn  implantation 

dose  of  7x10^^  cm’^.  The  threshold  current  vs.  the  dose  is  shown  in  Figure  5.3.  The  dots 
are  experimental  data  and  the  solid  line  is  the  dependence  calculated  with  the  same  disor- 
dering parameters  in  Chapter  3.  They  agree  well  verifying  our  prediction  that  disordering 
of  the  MQWs  by  Zn  greatly  reduces  the  threshold  current  of  the  SQW  laser. 

5.3  Groove  Etching  Bv  RIE 

In  etching  the  groove  with  reactive  ion  etching  (RIE),  the  availability  of  a high 
quality  mask  is  important.  In  a conventional  photoresist  mask,  a relatively  thick  layer  of 
photoresist  is  needed  and  the  edges  of  the  patterns  are  generally  rounded  due  to  the  aniso- 
tropic development  of  the  photoresist  in  the  photolithographic  process.  Masks  with 
rounded  edges  will  degrade  easily  during  the  RIE  and  lead  to  rounded  groove  walls.  We 
use  a tri-level  mask  (Figure  5.4)  in  which  two  layers  of  photoresist  are  used  with  a thin 
layer  of  Ti  sandwiched  between  them.  A pattern  window  in  the  top  photoresist  layer  is 
opened  by  a photolithographic  process  and  then  the  thin  Ti  layer  in  the  window  is  etched 
with  CI2-RIE.  With  the  remaining  Ti  layer  serving  as  the  mask,  the  bottom  layer  of  the 


2.  Photolithography  - top  PR 
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Figure  5.4  Preparation  of  the  tri-level  mask. 
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photoresist  is  etched  by  O2  plasma  etching.  A nearly  vertical  wall  can  be  obtained  with  a 
thick  layer  of  photoresist.  Another  advantage  of  the  tri-level  mask  is  that  the  bottom  layer 
of  photoresist  can  be  baked  at  higher  temperatures  to  have  higher  resistance  to  RIE.  To 
etch  the  GaAs/AlGaAs  system,  it  is  necessary  to  have  high  base  vacuum  to  prevent  the 
AlGaAs  from  being  oxidized  during  RIE. 

To  test  the  etched  groove,  we  used  a conventional  SQW  laser  structure  as  shown  in 
Figures  5.5  and  5.6  (Note  the  absence  of  the  MQWs).We  have  fabricated  several  lasers 
with  an  etched  mirror  and  an  etched  groove,  respectively.  Figure  5.5  shows  the  measured 
P-I  curve  for  a laser  with  an  etched  facet  and  a cleaved  facet.  We  estimated  the  reflectance 
of  the  etched  facet  by  comparing  the  threshold  current  densities  of  the  etched-mirror  laser 
and  a control  laser  using  Relation  (4.6).  The  control  device  had  the  same  material  compo- 
sitions and  layer  thicknesses  but  had  two  cleaved  facets.  We  used  ct=8  cm‘^  and  Rq=0.3. 
The  reflectance  Rg  of  the  etched  facet  was  estimated  to  be  about  0.13.  The  measured  P-I 
curve  for  a laser  with  an  etched  groove  is  shown  in  Figure  5.6. 

5.4  Fabrication  of  the  Tapered  Coupler 


In  this  section,  we  present  the  fabrication  process  for  the  tapered  waveguide  cou- 
pler. The  gradual  disordering  of  the  multiple  quantum  wells  in  the  tapered  coupler  section 
is  realized  by  implantation  of  nitrogen  through  a tapered  Si02  mask  and  subsequent 
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Figure  5.5  Measured  P-I  curve  for  a conventional 
SQW  laser  with  an  RTF  etched  facet. 
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Figure  5.6  Measured  P-I  curve  for  a SQW  laser  with  an 
RJE  etched  groove. 
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Figure  5.7  Illustration  of  e-beam  evaporation  of  the  tapered 
Si02  mask. 
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annealing.  The  e-beam  evaporation  of  the  tapered  Si02  mask  is  illustrated  in  Figure  5.7.  A 
motor  driven  knife-edge  mask  covers  the  sample  before  the  evaporation  starts.  As  the 
evaporation  proceeds,  the  knife-edge  moves  in  the  direction  of  the  tapered  waveguide  cou- 
pler leaving  a linearly  tapered  Si02  layer  on  the  sample.  The  length  and  thickness  of  the 
Si02  is  controlled  by  adjusting  the  knife-edge  speed  and  the  evaporation  rate  of  Si02.  In 
our  experiments,  the  knife  edge  speed  is  20  |xm/sec  and  the  evaporation  rate  used  is  16  A/ 
sec.  For  a 3 mm  long  coupler,  we  have  a 0-2500  A thick  Si02  taper.  The  sample  with  the 
tapered  Si02  mask  is  implanted  with  Nitrogen  at  an  implantation  energy  of  125  keV  with 

a dosage  of  5x10^^  cm'^  and  then  annealed  at  700  °C  for  40  minutes.  Figure  5.8  shows  the 
photoluminescence  spectra  of  the  sample  along  the  taper.  The  photoluminescence  peak 
due  to  the  MQWs  shifts  along  the  taper  indicating  the  MQWs  are  gradually  disordered 
along  the  taper.  Various  impurities,  including  Zinc,  Fluorine,  Oxygen,  and  Nitrogen,  have 
been  used  to  disorder  the  MQWs  in  the  tapered  coupler  section  in  our  experiments.  Nitro- 
gen was  found  to  provide  the  lowest  optical  loss.  Losses  due  to  the  nitrogen  doping  was 
estimated  by  measuring  the  losses  of  a nitrogen  implanted  and  annealed  MQW  waveguide 
without  the  Si02  cap.  The  measured  loss  was  ~6  dB/cm.  In  the  tapered  coupler,  the  losses 
due  to  nitrogen  varied  along  the  taper  from  the  maximum  (~6  dB/cm)  to  the  minimum  (~0 
dB/cm).  The  average  was  ~3  dB/cm.  This  is  significantly  lower  than  the  12  dB/cm  losses 
our  group  previously  obtained  for  the  Zn  disordered  tapered  interconnect.  Figure  5.9 
shows  the  measured  optical  loss  due  to  nitrogen  doping. 
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Figure  5.8  Measured  photolunninscence  spectra  of  the  sample  along  the 
taper.  The  numbers  represent  different  positions  along  the  taper. 
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Figure  5.9  Meaured  optical  loss  of  the  tapered  coupler. 
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5.5  The  Integration 

The  fabrication  process  of  the  integration  of  a SQW  laser  and  a tapered  waveguide 
coupler  (Figure  2.7)  involves  the  material  growth,  the  MQW  disordering  in  the  laser  sec- 
tion, taper  coupler  fabrication,  the  etching  of  the  groove,  and  the  laser  fabrication.  The 
process  is  summarized  below. 

1.  Material  growth  by  MOCVD.  The  layer  composition  and  thickness  sequences 
are  shown  in  Figure  3.3. 

2.  Zn  implantation. 

1)  Clean  the  sample  in  TCE,  aceton,  methanol  and  DI  water 

2)  Use  copper  tape  as  mask  to  cover  the  taper  section  to  prevent  unwanted  Zn 
implantation  in  that  section. 

3)  Implant  Zn  at  an  energy  of  165  keV  with  a dosage  of  7x10^^  cm'^. 

4)  Clean  the  sample.  If  necessary,  use  the  oxide  plasma  etch  to  remove  any  res- 
idue from  the  tape. 

3.  Si02  tapered  mask 

E-beam  evaporate  Si02  using  the  motor  drive  knife-edge.  The  Si02  evapora- 
tion rate  and  the  knife-edge  speed  are  16  A/s  and  20  |i.m/s,  respectively.  The 
taper  thickness  varies  from  0 to  2500  A. 
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4.  Nitrogen  implantation 

1)  Cover  the  laser  section  with  copper  tape. 

2)  Implant  Nitrogen  at  energy  of  5x10^^  cm'^. 

3)  Clean  the  sample.  Use  oxygen  plasma  if  necessary. 

5.  Annealing 

1)  E-beam  evaporate  500  A thick  Si02  on  the  sample 

2)  Place  the  sample  upside  down  on  a silicon  wafer  and  anneal  at  700  °C  for  40 
minutes. 

3)  Remove  the  Si02  cap  layer  in  buffered  oxide  etchant  (BOE) 

6.  Groove  etching 

1)  Spin  photoresist  S 1400-33  at  5000  RPM  for  40  sec.  Bake  at  120  °C  for  30 
minutes. 

2)  E-beam  evaporate  300  A Ti. 

3)  Spin  photoresist  S 1400- 17  at  5000  RPM  for  40  sec.  Bake  at  90  °C  for  30 
minutes. 

4)  Expose  the  groove  patterns  for  2.5  sec. 

5)  Develop  in  MIF  319  for  40  sec.  Rinse  in  DI  water  and  blow  dry  with  N2. 
Hard  bake  at  120  °C  for  30  min. 

6)  CI2  reactive  ion  etching  for  Ti. 

7)  Oxygen  plasma  etch  for  the  lower  layer  photoresist. 
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8)  Etch  the  groove  on  the  GaAs/AlGaAs  substrate  by  CI2  RIE.  Pump  the  sys- 
tem down  to  2x10’^  torr  before  starting  etching  to  prevent  A1  from  being  oxi- 
dized. 

9)  Remove  the  remaining  photoresist. 

7.  Metallization  for  the  laser 

1)  Spin  photoresist  S 1400-23  at  5000  RPM  for  40  sec.  Bake  at  90  °C  for  30 
minutes. 

2)  Expose  the  50  pm  stripe  patterns  for  4 sec. 

3)  Develop  in  MIF  319  developer  for  40  sec,  rinse  in  DI  water  and  blow  dry. 

4)  E-beam  evaporate  Ti  (200  A)  and  Au  (KXX)  A). 

5)  Lift-off  in  aceton. 

6)  Lap  the  backside  until  the  sample  thickness  is  ~1(X)  pm. 

7)  Clean  the  sample  in  warm  TCE,  aceton  and  methanol,  rinse  in  DI  water  and 
blow  dry. 

8)  E-beam  evaporate  1000  A thick  AuGe  on  the  backside  of  the  sample. 

9)  Cleave  the  sample  for  measurements. 

Utilizing  the  techniques  we  have  developed,  the  integration  of  a GaAs/AlGaAs 
SQW  laser  and  a tapered  waveguide  coupler  has  been  accomplished.  Figure  5.10  shows 
the  I-V  curve  of  the  laser  diode  of  the  integrated  structure.  The  P-I  curve  of  the  laser  is 
shown  in  Figure  5.11.  We  also  compared  the  power  output  from  the  laser  side  and  from  the 
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taper  side  at  various  injection  current  levels.  The  output  from  the  taper  side  is  about  5%  of 
that  from  the  laser  side,  which  is  much  lower  than  that  in  the  ideal  case  we  calculated  in 
Section  4.4.  Scattering  due  to  the  roughness  of  the  etched  groove  facets  and  diffraction  in 
the  groove  may  have  caused  much  of  the  losses.  Further  work  is  necessary  to  reduce  the 
roughness  of  the  etched  groove  and  hence  increase  the  output.  Another  possibility  would 
be  to  incorporate  a DFB  stmcture  for  the  laser  section.  Unfortunately,  it  comes  at  the 
expense  of  material  regrowth. 
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Figure  5.10  Measured  I-V  curve  of  the  laser  diode  in  the 
intgrated  structure. 
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Figure  5.1 1 Measured  P-I  curve  of  the  laser  in  the  integrated  structure. 


CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 


In  this  dissertation,  we  studied  theoretically  and  experimentally  a novel  monolithic 
integration  of  a single  quantum  well  laser  and  a tapered  waveguide  coupler  in  the  GaAs/ 
AlGaAs  materials.  Compared  to  the  existing  laser-waveguide  integration  structures  using 
butt  coupling  or  directional  coupling,  this  integration  scheme  avoids  complicated  and 
costly  material  regrowth  without  sacrificing  the  design  and  fabrication  flexibilities. 

The  key  technologies  in  this  integration  are  the  disordering  of  the  multiple  quan- 
tum wells  (MQW)  in  both  the  laser  section  and  the  tapered  coupler  section,  and  etching  of 
a groove  between  the  laser  section  and  the  tapered  coupler  section.  The  single  quantum 
well  (SQW)  waveguide  for  the  laser  and  the  MQW  as  the  output  waveguide  for  modula- 
tion are  vertically  stacked.  The  disordering  of  the  MQWs  in  the  laser  section  is  necessary 
in  order  to  achieve  low  threshold  current  density  for  the  laser.  We  numerically  modeled  Zn 
induced  disordering  of  the  MQWs.  We  first  used  an  explicit  finite  difference  method  to 
solve  the  diffusion  equation  in  which  the  diffusion  coefficient  is  concentration  dependent. 
The  resulting  numerical  Zn  diffusion  profile  as  a function  of  time  is  used  to  calculate  the 
(Zn)  concentration  dependent  interdiffusion  coefficient  of  A1  and  Ga  in  the  quantum  well 
and  barrier  layers.  Although  some  researchers  have  used  the  Boltzmann— Matano  method 
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and  the  error  function  approximation  to  investigate  the  Zn  induced  interdiffusion,  they 
ignored  the  time  dependence  of  the  interdiffusion  coefficient  and  the  two  methods  are  not 
applicable  to  the  problem  we  are  dealing  with.  We  developed  a finite  difference  approach 
for  solving  the  interdiffusion  profile.  An  implicit  method  was  used  to  ensure  the  solution 
stability.  The  Zn  diffusion  in  the  MQWs  alters  the  compositional  profile  of  the  MQWs 
resulting  in  a modification  of  the  refractive  index  profile.  We  continued  the  process  to  cal- 
culate this  modified  refractive  index  profile.  We  assumed  the  refractive  index  of  the  quan- 
tum well  is  equivalent  to  that  of  the  AlGaAs  alloy  that  has  an  identical  bandgap  since  the 
exciton  effect  contribution  to  the  refractive  index  is  negligible  when  the  optical  energy  is 
far  below  the  bandgap.  We  again  used  the  finite  difference  method  to  solve  the  Schoe- 
dinger  equation  for  the  energy  levels  in  the  disordered  quantum  well  and  then  calculated 
the  corresponding  refractive  index  profile  in  the  MQW/SQW  laser  structure.  The  disorder- 
ing of  the  MQWs  influences  the  laser  performance  through  the  optical  confinement  in  the 
SQW  active  layer.  With  the  modified  refractive  index  profile,  we  calculated  the  variation 
of  the  optical  confinement  factor  with  the  disordering  parameter  using  the  multilayer  stack 
matrix  theory.  Finally,  we  performed  the  gain  calculations  for  the  SQW  laser  and  investi- 
gated the  dependence  of  the  threshold  current  density  on  the  disordering  of  the  MQWs. 
We  predicted  a sharp  reduction  in  the  threshold  current  density  when  the  MQWs  are  disor- 
dered, which  agreed  well  with  our  experimental  results  obtained  later. 

The  etched  groove  between  the  laser  section  and  the  taper  coupler  functions  as  a 
partially  transmissive  mirror  providing  the  needed  feedback  for  the  laser.  The  reflectance 
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of  the  groove  was  investigated  by  using  a multilayer  matrix  method.  The  laser,  groove,  and 
tapered  coupler  were  considered  as  “layers”  in  the  light  propagation  direction,  with  the 
taper  being  further  divided  into  many  “layers”  so  the  taper  can  be  approximated  by  a series 
of  small  steps.  The  mode  index  of  the  fundamental  normal  mode  is  used  as  the  refractive 
index  in  each  layer.  The  reflectance  of  the  groove  was  thus  calculated  as  a function  of  the 
groove  width.  We  found  the  reflectance  oscillates  with  respect  to  the  groove  width  and  its 
value  can  be  significantly  larger  than  that  of  a cleaved  facet.  We  used  a simple  approach  to 
calculate  the  threshold  current  of  the  laser  with  an  etched  groove,  in  which  the  threshold 
current  of  the  etched  groove  laser  is  expressed  in  terms  of  known  parameters  of  a laser 
from  the  same  wafer  but  with  two  cleaved  facets.  The  threshold  current  can  be  reduced 
with  appropriate  groove  widths.  However,  the  dependence  of  the  threshold  current  on  the 
groove  width  was  rather  difficult  to  measure  due  to  the  nonuniformity  of  the  etched  groove 
and  the  scattering  associated  with  it.  As  such,  no  correlation  could  be  established. 

We  developed  key  fabrication  techniques  for  the  integration.  Zn  induced  disorder- 
ing was  realized  by  Zn  implantation  and  subsequent  annealing.  Implantation  was  chosen 
over  vapor  or  film  (deposited)  diffusion  because  of  its  ease  in  controlling  of  the  dose  in  as 
wide  a range  as  possible.  We  used  various  implantation  dosage  and  successfully  reduced 
the  threshold  current  from  over  1000  mA  to  -400  mA  by  increasing  the  implantation  dos- 
age indicating  the  disordering  of  the  MQWs  occurring  as  predicted  by  our  numerical  sim- 
ulation. The  groove  between  the  laser  section  and  the  taper  was  etched  by  CI2  reactive 
etching  (RIE).  We  used  a tri-level  mask  to  improve  the  mask  quality  and  ensure  vertical 
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walls  in  the  groove.  Regular  SQW  lasers  with  etched  facet  and  etched  groove  were  suc- 
cessfully fabricated  but  showed  nearly  twice  high  threshold  current  as  the  lasers  with  two 
cleaved  facets,  presumably  because  of  the  roughness  of  the  etched  facets.  Using  these 
techniques,  we  have  successfully  fabricated  a laser-taper  coupler  integrated  device 
although  the  threshold  current  and  the  loss  are  somewhat  higher. 

In  conclusion,  we  have  investigated  theoretically  and  experimentally  the  key  tech- 
niques for  the  integration  including  the  Zn  induced  disordering  of  the  MQWs  in  the  laser 
section  and  the  etched  groove,  and  demonstrated  the  integration  of  a SQW  laser  and 
tapered  coupler  in  the  GaAs/AlGaAs  material  system  without  any  regrowth.  Coupled  with 
other  research  efforts  in  the  Photonics  Research  Laboratory  on  the  tapered  waveguide  cou- 
pler as  well  as  in  electroabsorption  modulator,  this  dissertation,  in  effect,  completes  the 
integration  of  the  SQW  laser  - MQW  absorption  modulator  via  a tapered  waveguide  cou- 
pler in  GaAs/AlGaAs  without  any  regrowth.  Further  research  is  necessary  in  order  to  opti- 
mize the  groove  etching  process  as  well  as  improve  the  device  performance. 
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